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._ Résumé. — La méthode de dilution isotopique appliquée à la mesure de la section efficace de la 
réaction connue °!V{n, y)°?V, a donné des résultats très satisfaisants. Elle a donc été utilisée 
pour létude de la formation de Li dans 160. La cible irradiée par des protons de 155 MeV à Orsay 

Simultanément, la section efficace de formation de *Be dans les mêmes conditions a été déter- 
Le rapport de ces deux sections efficaces est en accord avec un calcul publié récemment par 


Abstract. — The isotope dilution method applied to the determination of the known cross 
section of the reaction °!V{n, y)°?V has given satisfactory results. 


It has then been used for 


The target, bombarded with 155 MeV protons is tri-distilled water. 
Simultaneously the production cross section of Be under the same conditions Was measured 


by counting the 480 keV y rays. 


The ratio of these two cross sections is in agreement with recent calculations made by I. Dos- 


trovsky et al. 


I. Introduction. —— Parmi les méthodes de 
mesure de la section efficace d’une réaction nu- 


* cléaire, celle qui consiste à déterminer le nombre 


de noyaux produits a surtout été développée par 
radiochimie. De ce fait, les réactions nucléaires 
conduisant à des isotopes stables ont été, jusqu’à 
présent, ignorées, à de rares exceptions près [1]. 


- Cette remarque, valable également pour les 1s0- 
topes à vie très longue, concerne notamment les 


réactions produisant ce qu’on appelle actuellement 
les fragments lourds (Li Be, etc.). Si la mesure de 
la section efficace de production de *Be dont la 


_ période et le rayonnement sont très commodes, est 


déjà effectuée dans certains cas [2], il n’en va pas 
de même en ce qui concerne Li et *Li. 

Nous nous sommes donc proposés d’entreprendre 
de telles mesures et, à cet effet, nous avons cons- 
truit un spectromètre de masse qui permet, à 
l’aide de la méthode de dilution isotopique [3], 
d'obtenir le résultat désiré. 


II. Le spectromètre de masse et la source 
thermo-ionique. — Ce spectromètre de masse [4], 


dont la figure 1 donne une vue d’ensemble, est un 
appareil à simple déviation magnétique dont le 
rayon est de 22 cm et le secteur magnétique de 600, 
Afin d'augmenter la sensibilité de la détection, le 
collecteur d’ions classique (cylindre de Faraday) 
a été remplacé par un multiplicateur d'électrons à 
dynodes de cuivre-beryllium, qui fonctionne avec 
un gain de 10 000 [5]. La figure 2 montre un 
spectre de masse du chrome naturel obtenu avec 
cet appareil. 

La source d’ions que nous utilisons est du typ2 
thermo-ionique. En effet, seule cette source est 
capable de produire des ions avec un rendement 
suffisant à partir de quantités d’éléments solides ne 
dépassant pas, dans certains cas, 10712 g et même 
105 

Elle nous a permis d’obtenir, à partir de 5.1077g 
de chrome se présentant sous forme d’acide chro- 
mique, un courant d’ions $?2Cr+ de 101? À pendant 
45 minutes. On obtient de même 1071? A pendant 
20 minutes avec 10719 g de lithium sous forme de 
sulfate. 

Cette source convient donc particulièrement bien 
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à l’étude des réactions nucléaires et spécialement à 
celles que nous avons envisagées et qui conduisent 


à des quantités de lithium de 1071 g à 101 g 
après un temps d'irradiation raisonnable. 


III. Mesure de la section efficace de la réaction 
51V (n, y) 5?V par la méthode de dilution isotopique. 
— Afin de valider la technique que nous allions 
utiliser pour doser le lithium, et obtenir une indi- 
cation sur les possibilités de cette méthode, nous 
avons mesuré la section efficace d’une réaction déjà 
connue [6]. 


PROCESSUS EXPÉRIMENTAL. — 65 mg de vana- 
dium naturel ont été irradiés pendant une semaine 
par un flux de neutrons thermiques de l’ordre de 
1,6.101? neutrons/sec/em? dans la pile P, à Saclay. 

La réaction 5V{(n, Y)52V conduit à la formation 
d’un isotope stable du chrome par émission $— : 


PEN Tire — 52Cr stable. 
3,7 min 


Pour effectuer le dosage du chrome 52, nous 
avons utilisé comme entraîneur du chrome 53 très 
pur obtenu à l’aide du séparateur d’isotope à 
double déviation magnétique d'Orsay [7]. La pro- 
portion de 5?Cr dans cet entraîneur était ramenée 
à 1/900. 


N°11 


DÉTERMINATION PAR SPECTROMÊTRIE DE MASSE 


TABLEAU I 


PROPORTIONS EN ATOMES (%) 


Isotopes 50 52 53 
Cr naturel 


Cr entraîneur 


(Orsay) 7 0,11 99,89 
Réactifs + Vanadium 

+ entraîneur 0,12 2,37 98,44 
Réactifs + Vanadium 

RE 0,11 3,68 96,15 


irradié + entraîneur 


Le chrome naturel introduit par contamination 
des réactifs chimiques et de la cible de vanadium 
donnait une contribution en 5?Cr qui a également 
été déterminée par dilution isotopique. Cette conta- 
mination correspondait à une quantité de 5?2Cr 
de 0,52 ug pour environ 30 g de réactifs et 65 mg 
de vanadium. 

Le traitement chimique qui précède l’analyse 
isotopique réduisant l’émission thermo-ionique 
chaque dosage a nécessité l’emploid’environ 20 pg, 
d'entraîneur, au lieu du microgramme qui est 
amplement suffisant quand on utilise de l’acide 
chromique pur. 

Le résultat de nos mesures au spectromètre de 
masse (résumé dans le tableau 1) a donné pour la 
quantité de ‘?Cr produite par l'irradiation la 
valeur : 0,30 0,02 ug. La valeur de la section 
efficace déduite des mesures ci-dessus est donc de : 


5,4 + 0,7 barns. 


L’incertitude relativement importante pro- 
vient pour moitié de la détermination du flux de 
neutrons ; d'autre part, la mauvaise émission 
thermo-ionique déjà mentionnée plus haut a con- 
tribué à l’augmenter. La valeur admise pour la 
section efficace est de : 


5,1 + 0,2 barns. (6) 


La méthode de dilution isotopique est donc utili- 
sable pour le dosage des isotopes stables produits 
par transmutation, mais il apparaît au cours de la 
mise en application de la méthode, certaines diffi- 
cultés que nous pouvons résumer à : 10 la conta- 
mination, qui limite véritablement la méthode ; 
20 l’extraction chimique de l’élément à doser qui 
conduit souvent à une mauvaise émission thermo- 
ionique, ce qui nécessite l'emploi d’une grande 
quantité d’entraîneur et, par suite, impose à cet 
entraîneur d’être isotopiquement très pur. 

Nous avons donc tenu compte de ces éléments 
dans l’étude suivante sur la section efficace de 


RaPPorT Poips Porps 
2 DE 52Cr D’ENTRAINEUR 
#Cr ug ug 
54 
2,38 
2,41 + 0,05 
0,07 cree 0,52 + 0,01 21,90 + 0,20 
É 100 
3,81 + 0,03 
0,06 ere 0,82 + 0,01 21,90 + 0,20 


production du lithium, en réduisant le nombre des 
manipulations indispensables et en utilisant une 
méthode ne nécessitant pas d'opération chimique. 


IV. Production de Li par la réaction 
Pisomev + 10. -—— La formation de fragnents 
lourds par bombardement de protons est actuel- 
lement interprétée de diverses manières. Dans le 
cas des noyaux lourds, il peut s’agir d’évaporation 
d’agglomérats (Be, Li, B, etc.) ou d’interaction 
directe, le proton incident cédant de l’énergie à une 
partie du noyau agglomérée pendant un court 
instant par suite des fluctuations de la matière 
nucléaire [9]. 

Dans le cas des noyaux légers, ces deux pro- 
cessus peuvent entrer en compétition avec la spal- 
lation du noyau bombardé, le noyau résiduel étant 
alors un fragment lourd. 

À la suite d’une suggestion de M. Radvanyi, nous 
nous sommes proposés d'étudier les variations en 
fonction de l’énergie des protons de la section effi- 
cace de production de fi, Là, Be dans 160 afin 
de déterminer la contribution de la spallation et 
celle de la formation directe. 

Nos premières mesures ont porté uniquement 
sur fL1 et Be. 

Le dosage de 7Be a été effectué par comptage 
des y de 480 keV. 

En ce qui concerne 6Li, c’est la méthode de dilu- 
tion isotopique que nous avons utilisée. La cible 
irradiée était constituée par un cylindre rempli 
d’eau placé à la sortie du synchro-cyclotron d'Orsay, 
les protons ayant une énergie de 155 MeV. 4 

Le choix de cette cible nous a permis d’utiliser 
une méthode extrêmement simple pour la récupé- 
ration des produits de réaction : aucune extraction 
chimique n’est nécessaire : l’eau irradiée est 
simplement concentrée jusqu’à une goutte qui est 
déposée sur le filament de la source d’ions, pour 
l'analyse isotopique du lithium. 

L’eau employée est tridistillée ; la quantité de 
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lithium résiduel, mesurée très précisément par dilu- 
tion isotopique (1). 

Lors des mesures, la quantité de lithium naturel 
contenue dans l’eau était de 107° g ; celle de 6Li 
formée après 10 heures d'irradiation d’un cylindre 
profond de 6 em était d’environ 1071 g. Le flux de 
protons, étalonné à l’aide d’une feuille d'aluminium 
était de 0 2uAA. 


Résultat. — Compte tenu du fait que du li- 
thium 7 est formé en même temps que Li, la 
ction efficace trouvée pour la formation de $Li 


dans 160 est de 
12 + 4 mb. 


Dans les mêmes conditions (énergie des protons 
incidents : 155 MeV ; perte d’énergie dans l’eau : 
30 MeV) celle de la formation de 7Be est de : 


5 + 1,5 mb. 


(2) Avec du lithium séparé à Harwell : ($Li) /(7Li) = 60. 


Nos 


Nous n’avons obtenu jusqu’à présent qu’une 
seule valeur concernant $Li, mais 1l est intéressant 
de noter que le rapport 


section efficace $Li 
section efficace *Be 


-est sensiblement égal à celui que donnent les calculs 


sur l’évaporation publiés récemment par Dos- 
trovsky et al. [10] et effectués pour des protons 
de 940 MeV et des noyaux de cuivre. Par ailleurs, 
si l’on considère qu’en plus de l’évaporation, la 
spallation de 160 intervient pour une grande part 
dans la formation des fragments lourds, nos résul- 
tats concordent avec ceux de E. Baker et al. [2]. 

Il nous reste donc à les confirmer en améliorant 
la précision des mesures et à doser le lithium 7, ce 
qui exige une cible d’eau extrêmement pure. 

Je remercie M. Bernas pour les conseils et l’aide 
qu’il m’a constamment prodigués au cours de ce 


travail. 
Manuscrit reçu le 23 juillet 1960. 
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MODÉLE UNIFIÉ 
POUR LES PHÉNOMÈNES DE PHOTODÉSINTÉGRATION NUCLÉAIRE VERS 20 MeV 


Par GEorGEs MONSONEGO, 


Laboratoire de Physique Nucléaire, Faculté des Sciences, Orsay (Seine-et-Oise). 


Résumé. — Nous étudions la résonance géante à l’aide d’une formulation unifiée qui rend 
compte du caractère individuel et collectif du phénomène. Nous faisons deux transformations 
unitaires sur l’hamiltonien du noyau, en introduisant deux nouveaux opérateurs conjugués æ&et 7, 
avec une condition subsidiaire sur les vecteurs d’état puisque nous vons étendu le nombre de 
degrés de liberté. Le nouvel hamiltonien comporte une partie collective en termes de œ et #, un 
partie intrinséque, et des termes d’interaction. Cete formulation, fondée sur le succès du modèle 
phénoménologique de Bohr et Mottelson, semble être satisfaisante ; les calculs théoriques de 
l’énergie de résonance, de la largeur de raie, de la section efficace, sont en accord avec l’expérience. 
La condition supplémetaire joue un rôle important dans le calcul de la largeur de raie. Nous 
étudions aussi l’effet des forces d'échange neutron-proton ainsi que la distribution angulaire des 
particules émises. 


Abstract. — We formulate an unified model for the study of the giant resonance, which takes 
into account the individual and collective character of this phenomenon. We perform two 
unitary transformations on the Hamiltonian of the nucleus, introducing two new conjugate opera- 
tors &æ and %, with a subsidiary condition on the eigenfunctions since we have increased the number 
of degrees of freedom. Inthe new Hamiltonian, there are : a collective part in terms of &« and x 
an intrinsic part and interaction terms. This treatment based on the success of the Bohr and 
Mottelson phenomenological model seems to be successful : the theoretical calculations of the 
energy, width and cross section of the giant resonance, agree with the experimental data. The 
subsidiary condition plays an important rôle in the calculation of the width. We study 
also the effect of neutron-proton exchange forces and the angular distribution of the emitted 
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particles. 


I. Introduction. — Les différents modèles nu- 
cléaires qui ont été proposés pour l’absorption de 
photons d’énergie de l’ordre de 20 MeV par les 
noyaux se présentent en général sous deux formes 
en apparence antagonistes : modèles collectifs [1] 
et modèles à particules indépendantes [2]. Aucun 
de ces deux types de modèles ne donne une des- 
cription totale satisfaisante de la résonance géante. 
Les premiers ne peuvent rendre compte des phéno- 
mènes d'émission directe, tandis que les seconds 
donnent en général des énergies de résonance trop 
petites. Peu de tentatives [3] ont jusqu’à présent 
été faites pour unifier le modèle collectif et le 
modèle des couches. Nous nous proposons d’étu- 
dier ici d’un point de vue purement méthodo- 
logique (*) un modèle unifié en partant d’une for- 
mulation générale due à Jean et Touchard [4]. 


II. Formulation. — Nous supposerons comme 
dans le modèle phénoménologique de Goldhaber- 
Teller que l'absorption dipolaire est due à loscil- 
lation collective de l’ensemble des protons par 


(*) Les calculs détaillés concernant cette étude seront 
publiés dans la thèse de l’auteur actuellement en prépa- 
ration. Lee APE 


rapport à l’ensemble des neutrons. Les variables 
collectives seront : 


A 
Re D: x Ve (1) 
k=1 
ANT 
Es 1 D Pr lr 
EN 
avec 
NE 1/Z pour les neutrons : k = 1,2 ...Z 
ET} 11N pour les neutrons k = Z + 1, ... À 


où Z + N — À nombre de masse et r est la dis- 
tance entre le centre de gravité des protons et le 
centre de gravité des neutrons. 


{r, p] = ir. 
Soit 
S Ph y 
Or 


l’hamiltonien du noyau, V est le potentiel nucléaire 
supposé indépendant des vitesses ; alors 


A vi = SP avec M = mZN]A. (2) 


Introduisons à priori deux variables subsi- 
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diaires &œ et x canoniquement conjuguées (qui vont 
jouer le rôle des variables r et p) avec la condition 
de contrainte «D — O0 sur un vecteur d’état ® 
de À. En effectuant sur À les deux transfor- 
mations unitaires 

77) 


1 L 
U== exp (; mr) et US =, exN (—; 
Phamiltonien du problème devient (5) : 

R = %, + AH (3) 


_ A PER SU LE 
ie = (45) : sm  owmt-D.{p, (p, HT)’ 


— JCintrinsèque + collectif 


AH= Ip, Hjat (4) 


\ (P. D. = partie diagonale de) 


, est l’hamiltonien d'ordre zéro où les mouve- 
ments sont découplés. Le coefficient de «? est la 
partie diagonale du double commutateur dont on 
calculera la valeur moyenne sur l’état de base 
intrinsèque. L'influence du terme de recul p?/2M 
sur À est négligeable. Un calcul simple montre que 
cela revient à remplacer la masse m du nucléon 
par une masse réduite telle que 


m*[m = AJ(A —1) 


AH est un terme d'interaction couplant le mode 
collectif au mode intrinsèque. Les termes non écrits 
sont la partie non diagonale du double commu- 
tateur et la série des autres commutateurs triple, 
quadruple, etc. 


La condition de contrainte devient : 


VE 0 (5) 


où Ÿ est un vecteur d'état de Je. 
En prenant pour V une somme d’oscillateurs 
harmoniques : 


J prend la forme simple. 


2 

Fe — (x—}) + + Mw?a.r +5 Mo CE UE) 

La fréquence collective est alors égale à la fré- 
quence intrinsèque de chacun des nucléons ; nous 
retrouvons le résultat de Brink. A cause de (5) le 
terme de couplage Mo? a.r est inopérant ; la 
condition de contrainte assure donc ici le décou- 
plage des mouvements. Ceci montre l’importance 
de cette condition sur les termes de couplage. En 
général les fonctions d’essai que nous choisirons 
pour Ÿ ne satisferont pas (5) mais le résultat pré- 
cédent montre que tant que nous travaillerons sur 
Phamiltonien d’ordre zéro &, où les mouvements 
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sont découplés, nous pourrons ne pas en tenir 
compte. 


IIT. Énergie de résonance. — 1° POTENTIEL 
SANS FORCE D'ÉCHANGE. — Nous supposerons le 


potentiel nucléaire purement central V — 2, V(Tav) 


où rw = [ra — r,| et pour simplifier les calculs nous 
supposerons un noyau à symétrie sphérique. De 
plus dans l’état fondamental intrinsèque, on consi- 
dérera la densité nucléaire p — 3/4rkR* constante 
dans le noyau, nulle à l’extérieur et les effets de 


surface seront négligés. Le coefficient € des a? 


devient en développant le double commutateur et 
en passant à l’espace des moments : 


c = —() EL) E (&i — t;)? 
jh v(k) 3 < eikkri-r) > dk (6) 


v(k) est la transformée de Fourier de V(r;;). € n’est 
différent de zéro que si t, £ t; c’est-à-dire si & est un 
proton et 7 un neutron. La seule contribution à C 
est donc la partie neutron-proton des forces nu- 
cléaires ce qui traduit l’hypothèse de départ. 
<etk(ri-r;) > est la valeur moyenne de etk(ri—r;) 
calculée sur l’état fondamental intrinsèque. Les 
calculs ont été faits pour un puits carré, un puits 
gaussien et un puits de Yukawa avec R = r, AY 
tri 2E010RLE Ce 

Les paramètres de profondeur et portée sont 
ceux qui résultent des expériences de diffusion nu- 
cléon-nucléon (6). Le cœur dur a été négligé puisque 
selon les idées actuelles sur la matière nucléaire seul 
l'effet des forces à longue portée intervient dans 
les phénomènes collectifs. Enfin nous avons pris 
Z = N — A/J2 et nous trouvons que pour ces 
trois puits l’énergie de vibration collective est 


koc & 39 ATUS 


qui donne un bon accord avec l’expérience. 


20 POTENTIEL AVEC FORCE D'ÉCHANGE. — Nous 
nous limiterons à une interaction de la forme : 
Vab) = Vasf(i — x) + xPor] = Vabll — x) + Vo (7) 


æ étant la fraction de force d’échange et Pay l’opé- 
rateur qui échange les positions des particules & 
et b. L’hamiltonien du noyau 


EEE avec Ve= DV 
a<b 
se transforme par U, en 


H; = H + De Va + Z Vi +h 
<< <i 


# 
#© 
dl 
2: 


où 
e ZA 
hk; —= 2 Vi; = % 5) Vé(zi x) Ty 


1 ZAN\2 
5m (2x) E Véte — st + 


—. en n’écrivant que la composante z pour simplifier 


la notation, nous négligeons le deuxième terme 
de k, qui est nul en moyenne — h, est une série de 
terme général proportionnel à (A/ZN r}" donc 
rapidement convergente et nous pouvons nous 
limiter à une approximation quadratique. Avec ces 
approximations, la deuxième transformation uni- 
taire U, sur À, donne un hamiltonien % où 
variables collectives et intrinsèques ne sont plus 
séparées. Pour étudier l'influence des forces d’é- 
change sur les paramètres collectifs il est suffisant 
d’avoir l’hamiltonien collectif d’ordre zéro J£ qui 
est obtenu en prenant J en moyenne sur le vecteur 
de base intrinsèque ; nous obtenons alors : 


2 1 
+30 + Ci) & + Con GR + Com n°? (8) 


à 
I 


= -sm(2y) <i> | 
ù <em> » 
Car = . 5) she [r. 2, Vis — zj)?| 


Cora = _ CS) < [. [r 2 Vii(si — 25)*]] pa 


Dans ces expressions z est un proton, 7 un neu- 


_tron. Comme nous négligeons le spin des particules, 


il n’est pas nécessaire d’antisymétriser la fonction 
d’onde intrinséque et nous choisissons un produit 
d’ondes planes. Les calculs ont été faits avec un 
puits carré et nous trouvons : 


CRE 
| 1 + 0,8x 


RAC OC 


M* (10) 


et pour la fréquence de vibration : 


Les paramètres de la dynamique collective sont 
donc directement rattachés aux grandeurs carac- 
téristiques des nucléons et des interactions nu- 
cléaires. En particulier la fréquence des oscillations 
collectives est parfaitement corrélée aux para- 
mètres des forces collectives déduits des expé- 
riences de diffusion. De plus les forces d'échange 
n’apportent pas de modification sensible à l'énergie 
de résonance. Leur effet revient à remplacer la 
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masse Ï de l’oscillateur collectif par une masse 
réduite M* donnée par (10), effet qui est compensé 
par une variation sensiblement équivalente du 
coefficient de la forme de rappel. 


IV. Section efficace totale intégrée — L’opéra- 
rateur de transition dipolaire électrique normalisé 
à un flux incident de un photon par cm? par sec 
est : 


es 


Ds. = eV9r he FE 


en prenant Oz pour direction incidente des pho-: 


tons ;,u = LE letrs, = T 


at ty), Dy 5e 


transforme par U, et U, en tenant compte de (5) 
Ah eV Fe ZA a. 


L’absorption du photon correspond à l’excitation 
dipolaire du mode collectif l’état de base intrin- 
sèque ne changeant pas dans l’absorption en se 
limitant à l’hamiltonien d’ordre zéro. L’absorption 
fait donc passer de l’état dr, ... ra) of) à 
l’état d(r; ... ra) p,(æ) et en négligeant les forces 
d'échange la section efficace totale intégrée est : 


ae 2e DEN EEE y (EU 
OR 2 Ke DoDulDi 22 = 27 


Si l’on tient compte des forces d’échange à cause 
de (10) : 


of == 6y(1 + 0,8x). 


Cette augmentation due aux forces d'échange est 
la même que celle donnée par Levinger et Bethe à 
partir de calculs fondés sur l’utilisation de règles 
de somme [7]. 


V. Largeur de raie des noyaux magiques. — Les 
mouvements collectifs étant couplés aux mouve- 
ments intrinsèques par l'interaction AH après 
absorption du photon il se produit des transferts 
d'énergie de la structure collective à la structure 
intrinsèque. Ce sont ces interactions résiduelles qui 
sont responsables de la largeur de raie. Le mode 
collectif excité est donc un état quasi stationnaire 
dont l’énergie est absorbée par le mode intrin- 
sèque. Nous devons donc considérer l’hamiltonien 
total (3) et nous désignerons par [n' > l’état où le 
mode collectif est excité le mode intrinsèque restant 
dans son état de base ; Æ,, est la valeur propre 
correspondante ; |n > est un état intrinsèque excité, 
E, la valeur propre correspondante [p > et [m > 
désigneront des états quelconques de %Æ,. Nous sup- 
poserons en outre que l'interaction AJZ est établie 
adiabatiquement depuis le temps £ = — co où elle 
est nulle jusqu’au temps { — 0 où elle vaut AA et 
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que de même elle décroît adiabatiquement jusqu’au 
temps { — + co où elle s’annule. Dans la théorie 
des perturbations dépendantes du temps que nous 
utilisons [8] cela revient à remplacer les distri- 
butions qui figureront dans les résultats par les 
fonctions dont elles sont les limites. Nous cherchons 
une solution de la forme : 


avec 
j l 
Hit) — exp G et) AH exp (-; it) 
Les solutions P, vérifient : 


à 
D AH D jeter 0 
m 


+ 1% Opn’ Ô(t) pour — © < { < + co (13) 
avec les conditions initiales : 

bel Re 0e 
bplt) = 0 pourt << Oet pour tout p bp(0) = 0 pour p # n°. 


En passant aux transformées de Fourier : 


1 +00 2 (Ep—EN 
Bt = —5 [ UE < Fa > A7? 
1 +co î /— Et ne 


les équations (11) deviennent alors : 
(E — Er) < plFsln > 
= 2 < plAHm> < mFel# > + Er (11/ 


et en posant : 
< plFein > = < plUrin > 
€ n'|Frln > Ê(E — Es) pour p Zn (13) 
où 
U(x) = — im dx),  at(æ) — 1 


nous trouvons en calculant < n'|Fyln' > par (11!) 
et en utilisant (12) ; 


0 
ble ren 2 (14) 


= ee. n'\AHIn > C(E—E,) < n\Ugln >. 
Nous trouvons aussi que 


Ug = AH + AHE(E — &, Py) Ur + is (16) 


où P,, est un opérateur de projection (P,;|n° > = Ü) | 

Le facteur périodique de (14) montre que A’E/2 | 
est une modification du niveau collectif excité. Le | 
facteur exponentiel décroissant indique que l/2 est | 
une mesure de la largeur de raie. On peut alors | 
calculer L'/2 et A'E en séparant les parties réelles | 
et imaginaires de (15) et si on se limite à l’approxi- 
mation AA pour Uz 


V AH > PT /2 | 
P s KrA4m> LE da 
É AE (En — En)? QE 
Ag SP ANS IEEE (18) | 


2 n+n En Er 


La théorie des perturbations indépendante du 
temps au deuxième ordre permet d’ailleurs de | 
retrouver ces résultats à condition de rendre compte | 
de la dissipation de l’excitation collective dans les | 
modes intrinsèques èn remplaçant le propagateur | 


JE —R par Î[E, — Ro Il en résulte 
que le deuxième moment nucléaire 
M = DE MANMERS 
nEn’ 
est une limite supérieure de l2/4. Nous sommes | 


obligés de tenir compte ici de la condition (5). | 
Pour cela limitons-nous au premier terme de AJ, 


- [p, H]a et prenons pour fixer les idées : 


1 A A 
V = 9mo? 3, ré. viavee dv = D ACIZS ee 
k=1 k=1 

alors 

a LA 4 
AH Sr P,=mMmo° > L: a 

ñ|= 2 k=1 

; À 
+5[P,ve= Moa.r+=[p, vla. 


Le premier terme est inopérant à cause de (5) et 
il en résulte que nous pouvons travailler avec les 
mauvaises fonctions d’ondes (ne satisfaisant pas à 
la condition de contrainte) à condition de rem- 
placer AA par 


ed 1 
AH ==[p, v| & 


les valeurs numériques de C et D sont tirées des 
schémas de Nilsson [9] ; les calculs ont été faits à 
partir de M, sur 160, #Ca, %7r et on trouve les 
largeurs 2,8 MeV, 3 MeV, 4,8 MeV, ordres de gran- 
deur tout à fait comparables aux données empi- 
riques qui sont objectivement 4 ; 4,2 et 5 MeV à 
1 MeV près. En repérant les énergies par rapport à 
l’état de base intrinsèque, l'énergie perturbée de 
l’état excité est : 


V2 = AG + AE 
ÂËE étant donné par (18). 
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Les valeurs calculées de AE pour les trois noyaux 
précédents sont toutes inférieures au MeV. Il est 


_ possible en partant de (18) de mettre en moyenne 


AE sous une forme différente. En supposant que les 
énergies en dénominateur varient peu avec les va- 
riables angulaires, nous obtenons la forme 


AE = < OGC M0 > 
ab 


en moyenne sur l’état fondamental intrinsèque, 
G étant une constante. Tout se passe done comme 
si l'énergie de perturbation AE était due à une 
interaction dipolaire additionnelle. Cette conclu- 
sion est tout à fait analogue au résultat 


 d’Elliott [10] qui construit des mouvements collec- 


tifs de rotation en supposant chaque nucléon dans 
un puits d’oscillateur harmonique et en ajoutant 
une interaction quadripolaire de la forme 


Z r$ r? Pa (cos Gi;). 
2] 


Cette formulation unifiée permet donc de rendre 
compte de la largeur de raïe et il faut souligner le 
rôle particulièrement important joué par la condi- 
tion supplémentaire sans laquelle on est conduit à 
des valeurs trop grandes sans sens physique. C’est 
là probablement la raison de l’échec de tentatives 
effectuées pour rattacher la largeur de raie à un 
couplage entre mouvements coopératifs et intrin- 
sèques (Fujita [3]). 


VI. Largeur de raie des noyaux déformés. — 
Soit l', la largeur des noyaux magiques que nous 
supposons aussi celle des noyaux sphériques. En 
effet dans ces noyaux les nucléons en plus des 
couches magiques contribuent à la largeur mais en 
même temps ils interdisent certaines transitions 
aux nucléons situés dans la couche immédiatement 
inférieure. Pour un noyau déformé nous pouvons 
donc écrire : 


TT AR 


AT étant la constribution due à la déformation. 
En considérant une déformation statique l’équa- 
tion de la surface est : 


R=R|! +2 (3 cos? 6 —1)] 


et les surfaces équipotentielles vérifient : 


re)? 


top B—=1/4—X 


IPavec 


0 


en supposant la déformation à volume constant. 
Aux coordonnées x, y, z du noyau, sphérique vont 
correspondre les coordonnées 81/2? x, B1/2 7, (1/8) z 
du noyau déformé et : 


1 
xeS — nr? + ; M oè [8 + où) + = s] 


ni 
2M fe 
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d’où deux fréquences longitudinales et transver- 
sales or et © 


hoc = Roll — À) 
oct = hoc + (À/2)] 
et AT = Hes — œe). Il y aura dédoublement 
quand AT © co qui est le cas de 181Ta : quand 


AT < T°,/2 on n’observera qu’un élargissement de 
la raie de résonance. En exprimant À à partir du 
moment quadripolaire Q, calculé en supposant la 
densité de la matière nucléaire constante dans le 
noyau : 

o Q 


AT —-—. 
4 ZRKÿ 


Les résultats numériques que nous obtenons sont 
en bon accord avec l’expérience. Notre formalisme 
est ici analogue à celui d’'Okamoto [11]. 


VII. Émission de particules par effet direct. — 
Dans notre formalisme cette émission se fait en 
deux stades : Absorption du photon par le mode 
collectif ; décomposition dans le mode intrinsèque 
avec excitation d’un nucléon émis dans le conti- 
nuum. On peut encore utiliser la théorie des pertur- 
bations dépendantes du temps [0 > sera l’état 
initial de valeur propre Æ, (E, est l’énergie du 
photon incident si on mesure les énergies à partir 
de l’état de base du noyau) ; |n' > est l’état collec- 
tif excité ; [nr > l’état du nucléon excité. L’hamil- 
tonien du problème est : 


Ro TAH TH, 
où À, est l’interaction du photon avec le noyau. 
Y P ï 
———— ZLN 
H, =eV2r ET œ 


en normalisant à un flux incident de un photon 
par cm? par sec. Le calcul de b,(t) amplitude de 
probabilité de l’état [0 > permet en particulier de 
calculer la section efficace totale d’absorption du 
photon. 


LANPNCIPCE (T/2r) 
nc À m ho (Ey — hoc)? + (T?J4) 


sy) = 27° 


et nous retrouvons la section efficace totale intégrée 
en intégrant sur toutes les énergies Æ,. 

Nous appellerons rl mg les nombres quan- 
tiques du nucléon dans son état fondamental, n, {, m 
les nombres quantiques qui correspondent au nu- 
cléon excité. NV, R, À\ sont les nombres quantiques 
collectifs et 


= dur) Z < 1Rm AIM > Yim(8, ®)INR À > 


est une fonction propre de J,. bu(r) est la fonction 
d’onde radiale du nucléon. Nous considérons le 
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phénomène pour { + co c’est-à-dire physiquement 
pour des temps # tels que l#/2% > 1 où on est 
certain que le mode collectif excité s’est décomposé 
dans le mode intrinsèque. À ce moment À étant 
le rayon de la voie de sortie pour r > À on peut 


écrire la fonction d’onde 


à pt 


p> =X2S dir); 00; Mm>e A w. (19) 


S, est l'amplitude pour l’absorption du photon 
avec émission du nucléon dans la voie « l », k est 
le nombre d’onde du nueléon émis. Pour r < À 


îi 
Ent 


lb > = bft—> co) QufÆr)l; 00; Im>e (19’) 
A partir de (19) et (19°) 
1 

Sp = ———— bnlt co) WlKR 20 
PEUR) (t— ce) Qu(ÆR) (20) 

où Æ est le nombre d’onde interne du nucléon. 

1 

Ci =TniZ lS1l?. (21) 
Il n’y a que deux transitions possibles L, — 1, +1 
ou Lo — lo — s la distribution angulaire totale 


calculée à partir de (20) contient des termes d’inter- 
férences entre ces deux types de transition et est 
de la forme a + b sin? 0. Ces termes d’interférence 
permettent au rapport b/a d’être négatif ce qui va 
dans le sens des expériences de Spicer sur 160. 
Par ailleurs les sections efficaces calculées au voisi- 
nage de la résonance sont : 
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où t — 1/Z ou — 1/N suivant que la particule 
émise est un proton ou un neutron. 


Dei = [ Ver des dr 


a été tabulé par Wilkinson [2] 


2 
D = EU ny RE GER) 


est la largeur d'émission de la particule, v,le facteur |} 
de pénétration. Dans tous nos calculs nous avons | 


idéalisé les mouvements du nucléon dans le noyau 
en l’assimilant à celui d’une particule plongée dans 


un puits de potentiel. Une confrontation avec l’ex- 


périence a été faite sur ?7AI (y, n (12). Nous 


trouvons un bon accord qualitatfi 


Conelusion. — Ce type de formulation unifiée 
présente donc un intérêt méthodologique certain. 
Ce modèle n’introduit plus empiriquement les para- 
mètres de la dynamique collective mais les exprime 
à l’aide de grandeurs caractérisant les mouvements 


des nucléons et leurs interactions. L’emploi simul- | 
tané du modèle des couches considéré comme |/h 
«structure modèle intrinsèque » et du modèle collec- ! 


tif pur permet l’interprétation et le calcul simple 
de la largeur de raie. Il permet aussi de retrouver 
des résultats équivalents à ceux de Wilkinson dans 
le phénomènç d’émission directe. Enfin il convient 
de souligner le rôle très important joué par la 
condition supplémentaire. 


Si >h+1 Cette étude m’a été suggérée par M. le Pr Maurice 
On Z NN ANT io IST Jean. Les discussions que nous avons eues tout au 
PT ee ] Us El Di +1 ZM +1 long de son élaboration m’ont grandement aidé. Je 
(En — Enr)? + tiens à l’en remercier vivement. Je remercie égale- 
OL —1 (21) ment le Centre National de la Recherche Scienti- 
107 e2 JZN \?  Ty1 hoc 1e fique dont les allocations m’ont permis de mener 
= —— — ( ] DE ce travail à bonne fin. 
BNIC\A sente 2 +1 
(Er D En?) 2 ET : Te 
k Manuscrit reçu le 30 juin 1960. 
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SUR LA THÉORIE STATISTIQUE DES CHAMPS NON LINÉAIRES 


Par JAPETERPETSKY, 


Professeur à l’Université de Moscou (1). 


Résumé. — On introduit une distribution de probabilités canoniques et microcanoniques pour 
un champ non linéaire classique arbitraire. On trouve les expressions de la densité fonctionnelle 
de probabilité pour les lignes d’univers qui représentent les solutions du genre particule des équa- 
tions non linéaires du champ. On montre que la formulation de Feynmann de la mécanique quan- 
tique n’est qu'un cas particulier d’une distribution canonique générale avec une température 
imaginaire. Cette température est interprétée comme la température statistique d’un thermostat 
formé de particules ayant des masses imaginaires. La constante de Planck représente alors une 
grandeur proportionnelle à la température du thermostat « imaginaire ». 


Abstract. — The canonical and microcanonical distributions of probabilities are introduced for an 
arbitrary classical, non-linear field. The expressions of the functional density of probability of the 
lines of u 1iverse are found, which give the particle-like solutions, of the non-linear field equations. The 
Feynmans formulation of quantum mechanics is only a special case of a general canonical distri- 


J 2%, od 
Ê = fcfvs de: ne, 


bution, with an imaginary temperature. 


1 That temperature is interpreted as the statistical tempe- 
rature of a thermostat constituted of particles having imaginary masses. 


The Planck constant 


is then a factor proportional to the temperature of the ‘‘ imaginary ” thermostat. 


I. La distribution microcanonique. — Exami- 
nons l’ensemble des champs classiques Y,(&,, x:, 
La, T3), LECT Lys Ty Ts +. Palo Ti Ta 
z,) dont les équations du mouvement sont obte- 
nues à partir du principe de moindre action : 


dx] 
Re 


drodradredre Sal) 


… où £ est le lagrangien du système, 5, est la valeur 


extrêmale de l’action S. 

En examinant du point de vue statistique cet 
ensemble de champs il est indispensable d’associer 
à chaque champ possible une certaine probabilité, 
a priori, ou plus exactement, une densité de pro- 
babilité dans l’espace fonctionnel. 

La plus simple des distributions de probabilité 
qui élimine les champs qui ne vérifient pas les 
équations du mouvement et qui donnent des 
valeurs équiprobables à toutes les solutions pos- 
sibles des équations du mouvement est évidem- 
ment la distribution. 

1È (2 


W = AlSn) 5 [Sn —s} UE 

où S est une fonction de Dirac, W et dF/ox sont 
l’ensemble de tous les champs et de leurs dérivées, 
mis sous une forme abrégée, À est une constante 
de normalisation qui est en général une fonction- 
nelle car SA dépend de l’énergie totale et des diffé- 
rents invariants intégraux (par exemple la charge 
totale). 

Par analogie avec la mécanique statistique clas- 


(4) Actuellement en mission au Service des Théories 
Physiques de l’Institut Henri-Poincaré. 


sique on peut appeler la distribution (2) distri- 
bution microcanonique. 

Il est évident que les valeurs moyennes de toutes 
grandeurs physiques F{Y(x) } sont définies par 
des intégrales fonctionnelles (?) 


F= [ F{ 4 }W{ ve) } apte). (3) 


Nous représentons par X l’ensemble des quatre 
coordonnées de Minkowi:ki æo, %3, To, 3 et 
nous utiliserons des accolades uniquement pour les 
arguments des fonctionnelles. Pour les arguments 
des fonctions nous utiliserons seulement des paren- 
thèses ou des crochets. 

La probabilité pour qu’une certaine fonction du 
champ Q[W(x)] ait une forme donnée g(x) sera 
évidemment déterminée par l’intégrale 


W { ata) } = [8 { ale — Qtvtæ)) } W{ vla } ave), 
(4) 
où D{y(x) } est une fonctionnelle à définie par 
l’équation : 


F{nla) } = 3 {vte) — nte) } F{ v(a) } dytai. (5) 


Il est évident que l’on peut déterminer cette 
fonctionnelle sous la forme 


5 { yla) } = f [exp if QnE(x) y(x) de | dE(x). (6) 


(2) Nous utiliserons ici pour l'intégrale fonctionnelle les 
notations que nous avons introduites en 1945 dans notre 
thèse de Doctorat publiée sous une forme abrégée en 
1949 [1]. 


ir, 


II. La distribution canonique. — Si un certain 
ensemble de champs y(x) parmi ceux qui ont été 
examinés ci-dessus (ou bien une certaine classe des 
solutions des équations du mouvement) peut être 
rattaché à un thermostat c’est-à-dire si l’on peut 
ne pas considérer les grandeurs utilisées dans ce 
problème physique alors les autres champs p(X) 
qui nous intéressent dans ce problème peuvent être 
donnés par une distribution canonique : 


W = Bexp [bS{ o(x)}] (7) 


où b est une constante numérique, B est une 
constante de normalisation qui est une fonction 
de b et des paramètres figurant dans S. 

De même, qu’en mécanique statistique la distri- 
bution (7) peut être obtenue à partir de (2) en 
utilisant quelques hypothèses simplificatrices natu- 
relles. 

En effet, si l’ensemble des champs Ÿ peut être 
partagé en deux sous-ensembles @ et y, et si l’on 
ne se préoccupe pas de la seconde, la densité fonc- 
tionnelle de probabilité pour les champs d’après (2) 
doit être de la forme : 


W{ ole)} = A{ Sn} f 5{5m — $ { pla), z{2) }1 ax(æ). 

(8) 

Pour un système de champs interagissant fai- 
blement on peut poser de façon approchée : 


S { plz), x(x) } = S1f o(x) } + S2{ x(x) } (9) 
alors, d’après (8), il vient : 
W { (x) } = DISn — S,{ œ(x) }] 


où (y) est une fonction déterminée par les proba- 
bilités du «thermostat » et qui est 


D(y) — A [oty—s;{ 1(æ) }] dx(æ). 


La fonction ® est évidemment donnée, à une cons- 
tante multiphicative près, constante qui est définie 
à partir des conditions de normalisation 


[wi ole) } delx) = 1. 


La fonction ® peut être évidemment calculée à 
partir de la formule (11) en se donnant la forme 
de S,f x(x) !. Cependant sa forme est déterminée 
univoquement si en plus de l’hypothèse (9) on 
admet que cette fonction à un caractère universel 
et ne dépend pas de la forme des fonctionnelles 
d’action des autres champs interagissant avec le 
«thermostat ». Dans ce cas, deux champs n’intera- 
gissant pas entre eux mais faiblement liés avec 
le thermostat (3), doivent vérifier la condition : 


DI(S" + 15") — SL qifx) } — S3f ox) |] 
= DIS — Si ex) } 1. DES” — S,f px) )] (13) 
() La liaison ne doit pas être nulle mais suffisamment 


faible pour que les conditions (9) soient approximati- 
vément réalisées. - . Re 


(10) 


(11) 


(12) 
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car les systèmes w.,(x) et p.(x) peuvent être consiihe 
dérés comme statistiquement indépendants. Maïh 
la seule solution de l’équation fonctionnelle (10! ) 
est évidemment une fonction exponentielle (7). || 


III. Les paradoxes de la théorie statistique des 
nous intéresse p(X) obéit à des équations du mou 
vement linéaires c’est-à-dire si le lagrangien £ esilh 
une fonction quadratique de ces champs et de leursh 
dérivées, la distribution (7) nous conduit à un para 
doxe que l’on appelle « la catastrophe ultra} 
violette » et qui fait que l’énergie du champ, même 
pour un domaine fermé, tend vers l’infini à causé Ë 
de sa distribution uniforme suivant les ondes mono: 
chromatiques de toutes les fréquences. Si lonk 
prend la distribution (2) au lieu de la distri} 
bution (7), le paradoxe demeure en ne changean 
que sa forme. Dans'le cas d’une distribution micro} 
canonique (2) ayant une énergie finie, l’énergie d 
toute onde monochromatique tend vers zéro et ne 
peut être différente de zéro que si l’énergie totale 
est infinie (4). 


æ1e 
œil 


IV. Les possibilités de la théorie des champs non 
linéaires. — Il en est autrement dans le cas no 
linéaire lorsque le lagrangien est par exemple 1æ} 
somme de termes quadratiques et de termes d’ordref 
pair plus élevé. Dans ce cas, pour une certaine# 
classe de lagrangiens, parmi les solutions ayant und 
sens physique et ayant une énergie, une charge, eteJh 
finies, 1l existe des solutions du genre particules 
(cf. par ex. : [2]) présentant des concentrations def 
champs stables (paquets). 
Pour de tels lagrangiens non linéaires, on peu 
bien entendu, avoir aussi des solutions sous formef 
d’ondes monochromatiques de faible amplitude, 
car pour de faibles amplitudes les équations d’ondesk 
deviennent linéaires. Cependant, les solutions quil 
existent dans le cas linéaire sous forme de sommesk 
d'ondes monochromatiques de toutes les fréquences 
et d’amplitudes approximativement égales (5) sont 
interdites dans le cas non linéaire car de telles! 
solutions supposent des concentrations du champ 
en forme de fonctions 3, de brève durée d’exis- 
tence et d'amplitude maximale aussi grande quel 
l’on veut. Ces dernières solutions sont impossibles! 
car, dans notre problème non linéaire, les seulesk 
solutions possibles avec une énérgie bornée sont 
des solutions du genre 
bornées du champ [2]. 


{ 


| 


| 
| 


particule, aves des valeurs} 


(*) Il est bien connu que ces difficultés sont évitées! 
en seconde quantification. Le but de cet article est l'étude! 
des possibilités de la théorie statistique des champs clas-! 
siques. Ce n’est qu’à la fin de cet article que nous mon-| 
trerons la relation entre théorie classique et quan til 

(5) Ces solutions sont la cause des « catastrophes ultra-| 
violettes » car si l’on avait une limitation du spectre du! 
côté de courtes longueurs d’ondes il n’y aurait pas non 
plus de divergences. | SUUTS 


=> 


>" 
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Pèdire (1), tft), ..., tft), 


[= Si nous examinons le cas d’un champ que rien 
[une limite dans l’espace et si nous considérons son 
Ménergie et certains invariants intégraux finis, seules 
les solutions du genre particule ont un sens phy- 
lu sique. Toutes les autres solutions sont de la forme 
[ù de paquets d'ondes quasi linéaires qui s’étalent 
ln rapidement sur tout l’espace et dont l’amplitude 
mtend par conséquent vers zéro. 


Ainsi, pour les lagrangiens non linéaires admet- 
tant des solutions du genre particule si l’on impose 


[que l’énergie et certains invariants intégraux 
soient finis, toutes celles 
‘après le principe de moindre action seront prati- 
[Mquement concentrées au voisinage des lignes d’uni- 
vers représentant le mouvement des maxima du 
penre particule. Mais dans ce cas la fonctionnelle 
d'action du champ peut aussi être représentée 
‘comme la somme des intégrales sur l’ensemble des 


qui sont possibles 


lignes d’univers de certaines trajectoires fonction- 
nelles qui s’appellent aussi en mécanique, fonction- 


nelles d'action. C'est-à-dire que l’on peut poser : 


San Î L{æ(t), æ’(t)] dt, 


où x(t) représente toutes les lignes d’univers, c’est- 
. et chaque x;(t) 
représentent les trois coordonnées d’espace en fonc- 
tion du temps d’une ligne d’univers k. 


(14) 


V. La probabilité des trajectoires des solutions 
du genre particules. — Il est évident que pour de 


- tels champs non linéaires la densité des probabilités 


peu être représentée non comme une fonction- 


«nelle de tout le champ mais comme une fonction- 


nelle des lignes d’univers représentant le mouve- 


ment des maxima des solutions du genre particule. 
D’après (2) et (14) nous obtenons la distribution 
microcanonique 


W { æ(t) } = A S[Sm— S{aæ(t)}] (45) 
et d’après (4) et (11) la distribution canonique 
W { æ(t) — B exp [bS { æ(t) 1] (16) 


comme distribution du système dans le thermostat. 
La valeur moyenne de toute grandeur F{ x(t) À 
sera évidemment déterminée d’après la formule : 


F = E { alt) } W { æ(i) } dæ(s). (17) 

A l’aide d’une formule analogue on peut égale- 

ment calculer la probabilité d’une valeur d’une 

grandeur physique à un instant donné. C’est ainsi 

par exemple que la densité de probabilité d’une 

valeur des coordonnées x à l’instant £ est donnée 
par la formule 


Wlx, t) — [3e — 2 W£atr)}de(r). (18) 


On obtient facilement, à partir de l’expres- 
sion (16) la densité de probabilité d’une configu- 
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ration de V points matériels : 1, Qo, ..., an Par 
exemple. Pour cela il faut intégrer (16) suivant 
l’ensemble des lignes d’univers æ(t) passant par la 
configuration de points donnée en cet instant. 
Autrement dit il faut calculer intégrale fonction- 
nelle (18). 

À l’approximation non relativiste le lagrangien 
du système est égal à la différence entre l'énergie 
cinétique et potentielle du système. En portant un 
tel lagrangien dans (16) et en intégrant par rapport 
à tous les æ({x) sauf un x(t,) donné on obtient : 


W(qu 2 -.., en) = exp [rb(F — U)], (19) 


où U est l'énergie potentielle du système, + est 
l’intervalle de temps minimum qu’il faut introduire 
pour effectuer pratiquement l'intégration fonc- 
tionnelle et Ÿ est une constante déterminée à 
partir des conditions de normalisation. Si l’on pose 
tb = 1/®O où O — ZT c’est-à-dire la constante de 
Boltzman multipliée par la température, l’expres- 
sion (19) est équivalente à la distribution canonique 
de Gibbs. 


VI. La mécanique quantique en tant que théorie 
statistique des systèmes ayant une température 
imaginaire. — Le cas particulier où dans la distri- 
bution (16) la constante b a une valeur imaginaire, 
présente un intérêt particulier. Si l’on met une 
telle distribution sous la forme 


W | æ(t) — B exp [is æ(t) }] 


on verra facilement que les valeurs moyennes des 
grandeurs physiques calculées à partir de la for- 
mule (17) avec la distribution (20) coïncide avec 
les valeurs moyennes quantiques de ces grandeurs 
si À est la constante de Planck divisée par 27. 

En effet, dans le cas non relativist?, l’intégrale 
de l’expression (20) prise sur toutes les trajectoires 
d’espace-temps passant par ce point et provenant 
de points antérieurs est égale à la fonction d’onde 
Ÿ(q1, I -.., Ian, t) comme cela a été montré par 
Feynman [3]. Tandis que l’intégrale de (20) prise 
sur toutes les trajectoires issues du point considéré 
à tous les instants ultérieurs est égale à D*(q:, 
Y2s +.) Qan) t). Mais d’après nos règles (16), (17), 
(18), l'expression (20) doit être forcément intégrée 
suivant toutes les trajectoires provenant de 
t — — coet allant à t — + oo, c’est-à-dire que la 
probabilité de la configuration donnée est égale au 
produit WY%, c’est-à-dire à la probabilité de la 
mécanique quantique. 

De façon analogue, ainsi que l’a montré Ria- 
zanov [4] on peut calculer la densité de probabilité 
de l'impulsion et des autres grandeurs physiques en 
supposant que l’expression (20) est la probabilité 
des trajectoires données partant de { = — co et 
allant à { — + co. Il est important de noter que les 
grandeurs physiques # dont les probabilités et les 
valeurs moyennes sont calculées à partir de (17), 


(20) 
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(18), (20), doivent être considérées comme des 
grandeurs classiques. 

Il n’est pas nécessaire d’introduire les opérateurs 
et les autres notions spécifiques de la théorie quan- 
tique orthodoxe. Toute la théorie peut opérer à 
partir de notions purement classiques et obtenir 
tous les résultats de la mécanique quantique. 

De sorte que tous les résultats de la théorie quan- 
tique sont obtenus comme conséquence d’un: 
théorie statistique générale du champ non linéaire 
classique ayant des solutions du genre particule si 
l’on admet que ce champ se trouve en équilibre 
statistique avec un thermostat ayant une tempé- 
rature imaginaire (puisque la constante b peut être 
considérée comme une grandeur inversement pro- 
portionnelle à la température). 

La seule différence entre cette théorie statistique 
et la théorie orthodoxe consiste en ce que la for- 
mule (20) n’est valable que dans la mesure où nous 
partageons la trajectoire en parties grandes par 
rapport aux dimensions de la particule (c’est-à-dire 
par rapport aux dimensions du maximum des solu- 
tions du genre particules). Pour des intervalles plus 
petits la probabilité des trajectoires perd sa signi- 
fication car il faut utiliser la probabilité d’une dis- 
tribution de champ donnée à l’intérieur de la par- 
ticule. Dans la théorie orthodoxe on considère que 
les postulats de la mécanique quantique sont 
valables pour des intervalles arbitrairement petits. 
Et par suite on admet que la formule (20) est rigou- 
reuse même pour des intervalles bien plus petits 
que les dimensions de la particule. 

De sorte que la différence entre la mécanique 
quantique habituelle et la théorie statistique du 
champ non linéaire avec un thermostat « imagi- 
naire » ne peut être déterminée expérimentalement 
qu’en étudiant la structure interne des particules 
élémentaires. 


VII. Signification physique de la température 
imaginaire. — On peut se poser alors deux ques- 
tions importantes : 

1° quel est le sens physique d’une température 
imaginaire ? 

20° qu'est-ce qui assure la constance de cette 
température ? 

Étant donné que la première question est fonda- 
mentalement physique, nous allons l’examiner 
dans le cas le plus simple. Examinons un gaz de 
particules relativistes de masses identiques m qui 
soit en équilibre avec le thermostat. 

La distribution de Gibbs (canonique ou micro- 
canonique) aura pour un tel gaz la forme : 


W(P;,P;, EN) e(mÈ ex) 


K 


(21) 


où PZ est l’impulsion d’une particule, 


Ex = VP£ + m° 


No 


son énergie (on a pris les unités telles que c = 1), 
et «y — 1/9 est un coefficient inversement pro- 
portionnel à la température. La fonction ® sera 
exponentielle pour la distribution canonique et une 
fonction delta pour la distribution microcanonique. 
Dans un autre système de référence, l’argument 
de la fonction © sera évidemment 


«D ex —B>DPr (22) | BK 
K K | 
où B et « constituent le quadrivecteur température |hs 
inverse. Pour des particules habituelles dont le | 
quadrivecteur (Px, ex) est du genre temps, le quadri- 
vecteur température inverse doit être, lui ausst, du 
genre temps, afin que la distribution prenne la | 
forme habituelle de la distribution de Gibbs dans M 
le système au repos. On peut cependant imaginer | 
des particules dont le quadrivecteur impulsion est 
du genre espace. Leur masse au repos |h 


m — V/e? — P? 
sera donc imaginaire. D’après (22), il parait alors {[h 
raisonnable de leur attribu2r un quadrivecteur 
température inverse, qui sera lui aussi du genre 
espace et donc d’introduire une température 1ma- 
ginatre définie par 


( s 
Eu 160 Vince” 


Pour de tels systèmes, la fonction de distri- 
bution f, c’est-à-dire la densité moyenne de parti- 
cules dans l’espace des moments peut être repré- 
sentée comme une fonction de (8, P), c’est-à-dire | 
comme une fonction de la température imagi- 
naire ($). 

Le problème de savoir s’il est possible d'admettre 
des particules de masse imaginaire a déjà été posé 
dans notre précédent article [5]. Nous avons 
montré que la théorie de la relativité et le deuxième 
principe de la thermodynamique n’interdisent pas 
l’existence de telles particules, cependant ces par- 
ticules ne peuvent pas transporter d’entropie 
négative comme le font les particules de masse 
réelle, 

Ainsi la « température imaginaire » peut être 
introduite pour décrire un ensemble statistique de 
particules ayant des masses imaginaires. Mais, |h 
comme nous l’avons vu plus haut, le coefficient | 
de la distribution (16) est inversement propor- |k 


(5) Il convient de noter que la fonction de distribution | 
ainsi obtenue est anisotrope. On peut cependant trouver |h 
une distribution isotrope qui serait la somme de distri- [M 
butions à températures inverses égales en module mais ayant | 
toutes les orientations possibles. Si toutes les orientations |N 
sont présentes avec un poids égal, la distribution obtenue | 
ainsi sera une fonction symétrique de B,P, c’est-à-dire, du [M 
module de l'impulsion divisé par la température imaginaire. 
Il importe de remarquer que pour une telle distribution, les 
moyennes de l’énergie et de l'impulsion sont nulles bien que | 
leurs moyennes quadratiques ne le soient pas. 
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. tionnel à la température du thermostat. Si celle- 


ci est imaginaire, le cocfficient b du système lié 
à ce thermostat, doit être aussi imaginaire. Par 
suite, nous supposons que le thermostat est un SYS- 
tème de particules de masses imaginaires ou bien, en 
utilisant la terminologie de la théorie des champs, 
nous dirons qu’on rapporte au thermostat toutes 
les solutions du genre particules qui ont un quadri- 
vecteur d’impulsion du genre espace, nous aurons 
alors un coefficient b imaginaire dans la distri- 
bution (16). 

L'hypothèse d’un « thermostat imaginaire » 


M explique aussi le fait que la température imaginaire 


est constante. En effet, d’après [5] les particules de 
masse imaginaire ne peuvent pas transporter d’en- 
tropie. Par suite l’entropie d’un système réel inter- 
agissant avec un « thermostat imaginaire » ne peut 
varier quelle que soit l’énergie du système réel. De 
sorte que l'hypothèse que le monde qui nous entoure 
est rempli de particules de masses imaginaires ayant 
une température imaginaire constante est admis- 
sible du point de vue physique et ne présente pas 
de contradiction interne. 

Il faut s’arrêter encore sur la question de la signi- 
fication physique et mathématique d’une densité 
de probabilité fonctionnelle complexe (20). Cette 
question s’est déjà posée dans [4] et on a essayé de 
Péviter en remplaçant les exponentielles par un 
cosinus. Cependant de cette façon on n’arrive pas 
à éviter des probabilités négatives. 

Un problème analogue se pose dans la mécanique 
quantique habituelle quand on introduit une den- 
sité de probabilité non dans l’espace de configu- 
ration mais l’espace des phases (cf. par ex. [7]et[8]). 
Là aussi la densité de probabilité est en général 
complexe ou négative. 

Il est évident que la densité de probabilité fonc- 
tionnelle (20) ne doit être considérée que comme 
une probabilité formelle ouguast probabilité, c’est-à- 
dire comme une grandeur qui permet, à l’aide des 
formules habituelles de la théorie des probabilités 


(ici la formule (17)), de calculer les valeurs moyennes 


ou les densités de probabilité des grandeurs qui sont 
mesurées directement par l’expérience. 

Une telle quasi probabilité peut être complexe ou 
négative si les probabilités des grandeurs physiques 
mesurables ainsi obtenues sont réelles et positives. 


. Cette dernière condition est évidemment vérifiée 
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pour les probabilités des grandeurs mécaniques 
habituelles calculables à partir des formules (20), 
(17), (18). 


VIIT. La quantification en tant que conséquence 
de la statistique relativiste d’un champ non li- 
néaire. — Ainsi tous les effets quantiques peuvent 
être considérés comme le résultat de l’interaction 
d’un système classique ayant une masse propre 
réelle avec un thermostat «imaginaire » c’est-à-dire 
avec un système formé de particule ayant des 
masses imaginaires et qui se trouve en équilibre 
statistique. La constante de Planck a alors la signi- 
fication d’une grandeur proportionnelle à la tempé- 
rature d’un thermostat « imaginaire » multiplhiée 
par une constante qui a les dimensions d’un temps. 

Cette constante de temps est égale à une cons- 
tante de longueur divisée par c. La constante de 
longueur est évidemment égale à celle qui figure 
dansles termes non linéaires des équations de champ 
(par exemple « ! » dans l’équation d’Heisenberg [9]) 
car elle définit les dimensions des solutions du 
genre particule c’est-à-dire les dimensions linéaires 
des particules élémentaires. De sorte que le nombre 
de constantes absolues dans cette théorie n’est pas 
plus grand que celui qui figure dans la mécanique 
quantique car la température du thermostat imagi- 
naire peut être considérée comme une grandeur 
relative dépendant de l’état de notre parcelle d’uni- 
vers. Dans les autres théories de la structure des 
particules élémentaires (cf. [9] par exemple) la 
constante de longueur est introduite comme une 
constante absolue complémentaire à côté de la 
constante absolue de Planck. 

De ce point de vue toute la mécanique quantique 
n’est qu’une théorie statistique relativiste du champ 
non linéaire, en tenant compte des états qui corres- 
pondent aussi bien aux particules ayant une masse 
réelle qu'aux particules ayant une masse imagi- 
naire. Les idées déterministes de Einstein et de 
Louis de Broglie se trouvent ainsi mises sous une 
nouvelle forme. 

En conclusion je voudrais remercier M. Louis 
de Broglie pour l’intérêt qu’il a porté à ce travail 
ainsi que Jean Pierre Vigier et Georges Lochak 
pour des discussions amicales fructueuses. 


Manuscrit reçu le 1€7 juin 1960. 
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Résumé. — L'énergie perturbée d'ordre n peut être exprimée à l’aide des seules fonctions 
perturbées d’ordre inférieur ou égal à n {2 pour n pair, ou (n — 1)/2 pour n impair. On généralise 
ainsi un résultat obtenu par Dalgarno et Stewart [1]. 


Abstract. — The perturbation energy of nth order can be calculated when the perturbation 


functions are known up to the order n/2 for n even, or (n — 1)/2 for n odd. 


This is a gerera- 


lization of a former result obtained by Dalgarno and Stewart [1]. 


Dalgarno et Stewart ont montré [1] que, pour 
calculer l’énergie perturbée d’ordre 2n + 1 corres- 
pondant à un hamiltonien de la forme A9 + AH%, 
il suffit de connaître les fonctions perturbées jus- 
qu’à l’ordre n. Nous allons généraliser ce résultat à 
un hamiltonien de degré quelconque en À. 

Considérons un système quantique régi par 
l'équation de Schr.dinger 


Hb = EŸ (1) 
ou, en notations de Dirac 
H(ÿ > = El 7, (1°) 


système pour lequel l’hamiltonien total Æ a la 
forme 


A HO) EH) LS HIDE, Me HO 
(2) 


À est le paramètre de perturbation et 4% désigne 

l’hamiltonien perturbateur d’ordre n. 
Conformément à la méthode des pertur- 

bations [2], nous écrivons l’énergie E sous la forme 


FR RO EUR ACIER Na EUR RO RE Ho ET) 


et la fonction d’onde sous la forme 
® — (0 + AD + A2U2 + ,,, + AM +... (4) 
soit encore, pour le vecteur ket correspondant 
dV> => +A> +2 > +... + nm > +... 
(47) 


En portant (2), (3) et (4') dans (1’),et en égalant 
les termes de même degré en à, nous obtenons, si 
nous voulons déterminer l'énergie perturbée 
d'ordre n, E®, le système de n + 1 équations 
successives 

d 
Z (HM— Em); m> = 0(j—0,1,...,n) (5) 


Ar = 0 


s 


soit, sous forme plus explicite 
(HO — E0)0 > — 0 
(HO — FÔ)|1 > + (HD — EU)|0 > = 0 


j==1 
(H(0) 2 E0))|; > de DA) — EG m)) |m 0) 
m= 0 
(5 ; 7) 
(HO — EU) in — j > 
n—i 
+ Y, (HG) — Etm)|n — ; HD = (5; n— 7j) 
m=1 
(H0) — Et0)]n — 1 > 
n—1 
+ Y (Htm = Etm)|n 1 Ho NT) (5: —"1) 
m=1 


n—1 
(AO — En > + Y (HU — Etm))|n — m > 
mMm=1l 


+ (H®) — Em)]0 > = 0. (5;7) 


/ 


Habituellement [2, 31, on multiplie scalairement 
les deux membres de (5 ; n) par le bra < O[, soit 


< 0[HD — En — 1 > 
+... + <'0[Hm-D — Em-0h > 
+ < 0[Htn) — Etm|g > = 0. (6) 


. On adopte ensuite la convention de norma- 
hisation 


AU ER N 
d’où l’on déduit les conditions d’orthogonalité 


< 0j > =0 qulqusoit 7 > 0. (7). 


. N°11 


Dans ces conditions, on tire de (6) 
En) = < OH — 1 > + < O[Hn — 2 > 
+... + < [Ho >. (8) 
On obtient ainsi le résultat bien connu que £° 
ne dépend que des fonctions perturbées jusqu’à 
- l'ordre x — 1 seulement, Nous allons montrer qu’il 
. est possible d'exprimer E"® à l’aide des fonctions 
. perturbées jusqu'à l’ordre n[2 si n est pair, à 
1" l’ordre (n — 1)/2 si n est impair. 
Montrons d’abord que l’on peut calculer E" à 
l’aide des fonctions jusqu’à l’ordre » — 2. Multi- 
. phons scalairement l’équation conjuguée de (5 ; 1) 
par In — 1 >,:l vient 


< 1[H0 — En —1 > 
+ < OH) — En —1> —0. (9) 
Éliminons l'élément matriciel 
< 1140 —_ Ep 1 > 
entre (6) et (9), nous obtenons 
EC A]HO En —_1 > 


RER UN NON 


DRE ENS 2 


n 
+ Y < 0jHm) — Elm)]n — m > = 0. (10) 
m=2 


Multiplions maintenant scalairement (5 ; » — 1) 
par le bra < 1|, soit 


< 1]H0 — En — 1 > 


n—1 
+ D  <1[H00 — Etmn —1— m> — 0. (11) 
m=l 


En éliminant encore le premier terme entre (10) 
et (11), on obtient 


n—1 
D <1/H0m) — Elm)]n — 1 — m > 


mM:=1 


n 
+ Y < 0[H00 — Elmin — m > — 0 (12) 
m= 0 


d’où l’on tire facilement 


n—1 
Em = YO < 14H00 — E)|n — 1 — m > 
m=1 


L 
Es < OH — Elmi|n — m> + < 6[HM]0 > (13) 
m=2 


c’est-à-dire que £® ne dépend plus que des fonc- 
tions perturbées jusqu’à l’ordre n — 2 (pour n > 2). 

L'expression (12) ne contient plus que les fonc- 
tions perturbées jusqu’à l’ordre 7 —2 et four- 
Dit m»Æ% si, de manière générale, nous con- 
venons de désigner par w_#»£% l’énergie £® expri- 
mée à l’aide des fonctions perturbées jusqu’à 
l’ordre n — j seulement. Nous allons montrer que 
Jon a une relation générale de la forme 


31 n—k 
D  <kAM — Emn—k—m> = 0 


k=0 m=j—k 


(14 ; j) 
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qui ne contient plus que les fonctions jusqu’à 
l’ordre r — j et fournira donc l’expression 4-2. 
Remarquons que (6) est un cas particulier de (14 ; j) 
pour 7 — 1, de même que (12) pour j = 2. 

Raisonnons par récurrence et montrons que l’on 
passe bien de (14 ; j) à (14 ; j + 1) en changeant j 
en 7 + 1 dans (14 ; j). Cette relation peut encore 
être écrite 


j—1 
D on Re pe 


CRD NE), 


<< k|HO — Etm\n — k— m> = 0. 
k=0 M=ÿ—k+1 


(14 ; j) 
En multipliant scalairement l’équation conjuguée 
de (5 ; j) par le ketln — j >, nous obtenons 
FA 
D < AH — EG—b|n — j > 
k=0 
LH < HO — En—j>—=0 (15) 


En multipliant scalairement l’équation (5 ; n —j) 
par le bra < j|, il vient 


< jlAU) — En — 7; > 
cn s << j|HU — En — j— m> = 0. (16) 
m=1 
En ajoutant (14; j) et (16) puis en retran- 


chant (15), nous obtenons finalement 


less . . 
>, < jlAG0) — En — j — m> 4’: 7 +1) 
m=1 
j—1 n—T% 
+ Y D << kJHM — En —k— m> = 0 


k=0 m=i—k+1 


ou encore 


j n—k 
D 5  <HHM— Emn—k—m> = 0 
k=0 Mm=i—k+1 


(LES EEE 


qui est bien la formule obtenue à partir de (14 ; 7) 
par la transformation 7 — 7 + 1 ; elle ne contient 
plus que les fonctions d’ordre inférieur ou égal 
à n — j — 1. Puisque la formule (14 ; j) est valable 
pour j = 0, 1, elle est par suite valable pour toute 
valeur de j compatible avec la valeur de n. 

De (14 ; j), on déduit facilement la valeur 


j—1 n—k 
o—nE = NX << k|HmM — En —k— m > 
k 


=1mMm=j—k 


—1 
ce S < 0[H6) — Etm|n — m > + < 0[Hm]0 >. (17) 


Mm=) 


Si l’on utilise les conditions (7), les termes en £” 
de la seconde somme sont tous nuls. 
Nous pouvons continuer ainsi Jusqu'à épui- 


EE 
99 
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sement de toutes les paires d'équations (5 ; 7) et 
(5; n — j). En remarquant que l’on élimine la 
fonction d’ordre n — j + 1 entre les équations 
(6:j—1)et(5;n— 7j +1) pour aboutir à l’ex- 
pression #£‘, nous voyons que nous devons 
avoir finalement 

— pour n pair : N— fmax. + 1 = (n/2) + 1d’où 
1mas —= n [2 et 


(n/2)—1 n—k 


r12)E 0) Fe D 2 
k=1 m=(n/2)}—k 


< k|HO0 — Ein — km 


ML cs 
+ XD  < 0/40 — Et)]n — m > + < 0[4m0 > 


m=m/2 
(18-1) +<l1jne ll ax + <i—21H® 152) 
n +1 
— pour R impair: N — Jmax. + À = ti STORE 
de _ 2, D EM < kn—k—m > (20-1) 
REC n ny 
d’où mers 5 et paire < 
se — pour # impair Î 
EN, AK n —1 n—1 | 
7 AN UUE— >; >» << k]HG) — Etn)|n — RES re nn AUD - | HG : s 
(es) k=1 ne +1 y (9 2 2 
: À on Dre 
+ CS. <oHm—_Emhn—m> + < 0[HM]0 > : : 
me FRE re Pa 
2 2 
(18-2) n—1 
\ $ EE n—k 
Dans la plupart des applications, on a seule- à > EU << kln — k— m 7. (20-2) 
ment 0e 0 FANAIER OEM 0 pour: 222; gs =" 
Manuscrit reçu le 15 juillet 1960. 
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L'expression (17) devient alors (pour 7 > 2 etl} 

j > 2), compte tenu de (7) 

(n—nEM = < j—1A|HOn — j7 
LC AlH ne jte 
MENÉS IMORE Se 


j—1 n—k 


PARU Pà 


k=1 m=j—k 


Em < k]n— k— m7, 


tandis que les expressions (18) deviennent : 
— pour 7 pair 


DE = < 5 — 1 | HO) 3 > 


[2] Par exemple Kramers (H: A.), Quantum mechanics 
(North Holland Publ. Cc), 1956, p. 188 sq. 
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RÉALISATION ET ÉTUDE D’UNE CHAMBRE A FISSION A BISMUTH, 
DÉTECTEUR DE NEUTRONS DE GRANDE ÉNERGIE 


Par L. KOWALSKI (*), 


Laboratoire de Physique Nucléaire de la Faculté des Sciences, Orsay. 


Résumé. — On a réalisé une chambre à fission multiplaque à bismuth, détecteur à seuil pour 
les neutrons de plus de 50 MeV. Cette chambre est utilisée comme moniteur du faisceau de neutrons 


de 143 MeV du synchrocyclotron d'Orsay. 


Abstract. — We have constructed à multiplate bismuth fission chamber which constitutes 


a threshold detector for neutrons of more than 50 MeV. 
the neutron beam of the Orsay synchrocyclotron. 


1. Introduction. — On sait que les noyaux de Bi 
bombardés par des particules de très grande énergie 
peuvent subir la fission [1{, 2, 3]. Ce processus a été 
utilisé pour détecter des neutrons d’une énergie 
supérieure à 50 MeV [4, 5, 6, 7] à l’aide de chambres 
d’ionisation à impulsions comportant des électrodes 
recouvertes de bismuth. 

Les fragments de fission émis par le bismuth sont 
beaucoup plus iomsants que les autres particules 
chargées présentes dans les faisceaux et on peut les 
détecter séparément en utilisant un discriminateur 
électronique. 

L'avantage spécifique des chambres à fission à bis- 
muth réside dans le fait que ce sont des détecteurs à 
seuil élevé qui peuvent être employés en présence 
d’un fond important de rayonnements ionisants. 

Les chambres à fission sont en général utilisées 
comme moniteurs et placées dans le faisceau direct 


de neutrons ; elles ne sont sensibles qu’aux neutrons 
l École 


(*) Habituellement à Polytechnique de 


Varsovie. 


This chamber is used as a monitor of 


énergiques effectivement utilisés dans les expé- 
riences. Pour des intensités très différentes, les taux 
de comptage restent proportionnels aux flux de 
neutrons. C’est pourquoi nous avons réalisé une 
telle chambre qui est utilisée sur le faisceau de 
neutrons de 143 MeV du synchrocyclotron d'Orsay. 


2. Description de la chambre à fission. — Pour 
construire une chambre à bismuth, on part des 
mêmes principes généraux que pour les chambres à 
fission à uranium couramment utilisées pour détec- 
ter les neutrons lents [8, 9, 10, 11]. Cependant, la 
faible section efficace de fission du Bi et le grand 
nombre de réactions en compétition augmentent 
les difficultés pratiques dans le domaine des grandes 
énergies. En effet, la section de fission 64 du Bi 
est égale à 36 mb pour des neutrons de 120 MeV 
et atteint environ 50 mb à 150 MeV [12, 13]. Pour 
réaliser un détecteur suffisamment efficace, nous 
sommes donc obligés d’augmenter le nombre de 
plaques recouvertes de Bi. Cela augmente la capa- 


Fic. 4. — Vue de l’intérieur de la chambre à fission, 
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cité de la chambre, qui fonctionne sur le principe 
des chambres d’ionisation à impulsions, et le fond. 
Ce fond est le résultat de l’empilement des petites 
impulsions dues aux diverses particules chargées 
provenant des nombreuses réactions nucléaires 
dans les matériaux de la chambre et c’est 
pourquoi le fond dépend de l'intensité instantanée 
du faisceau et par conséquent de la structure en 
temps [21]. 
La figure 1 montre l’intérieur de la chambre. 


N° cl 


Les électrodes en Al d’une épaisseur de 0,5 mm: 
sont supportées par des tiges métalliques reliées à 
l’enceinte par l'intermédiaire d’isolateurs en plexi-; 
glass. L’enceinte en fer reliée à la masse Joue le) 
rôle d’anneau de garde. | 

Les électrodes portées à la haute tension (fig. 2)) 
sont couvertes par des couches de bismuth. Le? 
parcours moyen des fragments de fission dans Bi! 
est égal à 4,3 mg/em? [14] et nous avons choisi une? 
épaisseur pour chaque couche de 0,64 mg/em?. Les; 


560kQ 


F1G. 2. — Schéma de la chambre et du montage électronique. 


SNS 


\Cible (CH) 


* ve a pe 


* De 
\ Se 
Chambres à 
a Bi 
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N 


F1G. 3. — Disposition générale de l'expérience. 


couches étaient préparées par évaporation de Bi 
sous vide et ces épaisseurs ont été mesurées à l’aide 
des rayons X [15]. La chambre est remplie par un 
mélange À +2 % N, qui est relativement insen- 
sible aux petites impuretés [16]. Nous avons utilisé 
une pression d’une atmosphère pour avoir une 
bonne stabilité de fonctionnement. 

En tenant compte des parcours des fragments de 
fission dans l’argon et des courbes de ralentis- 
sement [17], nous avons choisi une distance entre 
les électrodes de 1,5 em. La chambre contient 0,61 # 
de Bi distribué sur une surface totale de 955 cm? 
(18 couches) et sa capacité est égale à 60 pF. Nous 
avons pris une surface utile pour chaque électrode 
(surface couverte de Bi) de 53 em? de manière à 
couvrir entièrement le faisceau de neutrons du 
synchrocyclotron tel qu’il est défini par les diffé- 
rents canaliseurs. On obtient facilement des taux 
de comptage de l’ordre de 1 coup par seconde (le 
mouvement propre état négligeable). L’ensemble 
de l’appareillage électronique utilisé lors de cette 
étudeest représentéschématiquement sur la figure 2. 

Le temps de montée des impulsions est 
de 0,6 us. Pour diminuer la possibilité d’empi- 
lement, on a utilisé une faible résistance de charge 
(56 kQ) et une mise en forme dans le préampli- 
ficateur. Les impulsions à la sortie de ce dernier 
ont une forme symétrique et durent environ 2 y sec. 


Eneus 


Pour analyser le spectre d’impulsions, une porte 
linéaire, ouverte seulement par les impulsions dé- 
passant un certain seuil, est utilisée avant l’analy- 
seur à 256 canaux de manière à éviter sa surcharge. 
La chambre, lorsqu'elle est utilisée comme moni- 
teur, fonctionne seulement avec un amplificateur et 
un discriminateur suivi d’une échelle. 

Une petite source de particules « se trouve cons- 
tamment à l’intérieur de la chambre ; elle peut être 
tournée vers les électrodes pour un contrôle per- 
manent du fonctionnement de l'appareil. 


3. Étude du fonetionnement de la chambre à 
fission. — La figure 3 représente la disposition 
générale utilisée. 

Les protons accélérés dans le synchrocyclotron 
jusqu’à une énergie de 157 MeV tombent sur une 
cible de beryllium d’une épaisseur de 10 mm. Des 
neutrons énergiques sont projetés en avant à la 
suite d’une réaction d’échange (p, n) dans cette 
cible. Le spectre des neutrons possède un maximum 
vers 443 MeV [18]. 

Après être passés par un double canaliseur rec- 
tangulaire 30 X 50 mm? percé dans deux cylindres 
de fer de 50 cm de longueur chacun sceilés dans 
des blocs de béton lourd du mur de protection, puis 
par une ouverture circulaire percée dans un bloc, 
les neutrons traversent une cible en polyéthylène 
et arrivent sur la chambre à fission. Pour l’étude 
de cette dernière, le télescope à scintillateurs réalisé 
par P. Radvanyi et J. Genin [19] a été utilisé 
comme moniteur du faisceau en comptant à 160 les 
protons de recul provenant de la cible de (CHo)n. 
La figure 4 montre la vue extérieure de la chambre 
à fission et du télescope. 

Les flux de neutrons de plus de 20 MeV ont été 
déterminés à l’aide de la réaction Cl?(n 2n)C11 [20] 
qui a servi à étalonner la chambre. 

Pendant ces mesures, le synchrocyclotron a été 
employé aussi bien dans des conditions normales 
(450 impulsions de faisceau par seconde) que dans 
le régime dit « stochastique », c’est-à-dire avec le 


RÉALISATION ET ÉTUDE D’UNE CHAMBRE A FISSION 781 


— Chambre à fission et télescope utilisé lors des mesures. 


programme auxiliaire d’accélération améliorant 
considérablement la structure en temps du fais- 
ceau [21]. 


4. Caractéristiques de la chambre, — La figure 5 
représente Îe taux de comptage en foncon de 1: 
haute tension appliquée. La chambre est utilisée 


monit 


HT. 
1500 V 


250 500 750 1000 1250 


Fic. 5. — Nombre d’impulsions 
en fonction de la tension appliquée à la chambre. 


couramment avec une tension de 800 V suffisante 
pour recueillir les charges produites dans le gaz par 
un fragment de fission (champ électrique de 
530 V/em). Une tension plus grande ne provoque 
pas une diminution appréciable du temps de montée 
des impulsions. 

La figure 6 représente les spectres de hauteur 
d’impulsion provenant de la chambre à bismuth 
pour trois régimes de faisceau différents. 

La chambre utilisée comme moniteur fonctionne 
seulement avec un amplificateur-diseriminateur et 
une échelle (fig. 2). La caractéristique du seuil de 
discriminateur (courbe intégrale) est tracée sur la 
figure 6 en pointillé. Nous travaillons aux environs du 
point À sur la figure où l’influence des instabilités 
de l'électronique sur le taux de comptage est mini- 
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mum. En augmentant l'épaisseur de bismuth sur 
les électrodes, la chambre devient plus sensible, 
mais la partie plate de la caractéristique du seuil 
disparaît complètement, ce qui montre que la plus 


onit 


\ 
, : (EE ne 
Cm sw 
15 20 25 30 35  L0 
(Seuil de discr.) 


FrG. 6. — Spectres des hauteurs d’impulsion de la chambre 
à Bi pour différents régimes de faisceau. En pointillé, 
nombre d’impulsions en fonction du seuil du discri- 
minateur. 


grande partie du fond visible est due aux réactions 
nucléaires provoquées dans le Bi. 

L’épaisseur des couches de 0,64 mg/em? est un 
compromis entre la stabilité du moniteur et sa sen- 
sibilité aux neutrons. La caractéristique du seuil 
montre qu’une instabilité de lélectronique 
de + 5 % (gain de l’amplificateur et seuil du dis- 
criminateur) donne une variation du taux de comp- 
tage de + 2 %. On peut facilement changer les 
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plaques de Bi pour adapter la chambre à des exi- 
gences différentes. 

Le rapport du taux de comptage de la chambre 
à Bi à celui du télescope servant de moniteur pour 
ces mesures reste constant dans le domaine de varia- 
tion du flux de neutrons fourni par le synchrocy- |: 
clotron. 

Nous avons vérifié que la chambre est insensible }h 
aussi bien aux neutrons de 14 MeV qu'aux rayons y |M 
du Co. 

Dans le cas des neutrons, la chambre a été placée 
à une distance d'environ 50 cm d’une cible de | 
tritium bombardée par des neutons (— 10°n/s.4x); 
dans le cas du cobalt, elle était irradiée par le |h 
faisceau de y d’une source de 100 curies placée à 
une distance de 50 cm (= 20 r/h). | 

Les caractéristiques de la chambre sont assez 
stables ; en particulier, il est possible de l'employer | 
pendant plusieurs mois sans changer le gaz de rem- |E 
plissage. Pour faire le vide préalablement dans la |} 
chambre, il suffit d’utiliser une pompe à palettes. 
Aucune purification supplémentaire du gaz prove- | 
nant de bouteilles commerciales n’est nécessaire 
pour cette utilisation de la chambre. 
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M. Reïde pour l’aide qu’il m’a apportée dans la 
partie électronique de ce travail. 

Ma reconnaissance émue va au Pr F. Joliot qui a 
bien voulu m'accueillir dans son laboratoire. Je 
tiens à remercier M. le Pr Teillac pour l'intérêt qu’il 
a constamment apporté à cette étude, ainsi que les 
autorités françaises et polonaises qui m'ont apporté | 
leur aide matérielle. 
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DIFFUSION DES NEUTRONS LENTS PAR L'HYDROGÈNE LIQUIDE 


Par G. SARMA, 
Service de Physique du Solide et de Résonance Magnétique, C. E. N., Saclay. 


Résumé. La section efficace différentielle de diffusion des neutrons lents est mise sous une 
forme simple qui fait intervenir d’une part un facteur de forme provenant des éléments de matrice 
de transition, d’autre part une transformée de Fourier d’une fonction de corrélation n’intéressant 
que les centres de masse des molécules et caractérisant la dynamique du liquide proprement dit. 
Il est montré que la partie qu’on peut appeler « cohérente » de la diffusion, c’est-à-dire celle qui 
est associée aux corrélations entre molécules différentes, est négligeable. La diffusion fait essen- 
tiellement intervenir la fonction d’autocorrélation. L'importance du recul est mise en évidence. 
On montre qu'un modèle de gaz parfait tenant compte du recul est raisonnable dans le cas général. 
Il conduit à des raies inélastiques bien séparées. Le cas de la diffusion élastique aux petits angles 
doit permettre la mesure d’un coefficient de self diffusion. “ 


Abstract. — The differential slow neutron scattering cross section by liquid hydrogen is expres- 
sed, under some reasonable assumptions, as a product of two factors of simple physical meaning. 
The first factor, throughout spin correlations, involves transitions between molecular states and 
finally leads to a kind of an angular dependent scattering length, which is computed exactly for 
every transition in the rigid rotator approximation. One is then left with a second factor depen- 
ding only on the translational motion of the molecules, which turns out to be a Fourier transform 
of the well-known pair correlation function for a monatomic liquid. 

It is shown that what one might call the ‘‘ coherent ” part of the scattering, ï.e. the part invol- 
ving correlation between pairs of distinct molecules, is negligible. 

So the scattering is mainly related to the so-called self correlation function. Recoil appears 
naturally as a prominent feature. It is shown that the scattering is generally well described by 
means of a perfect gas model, taking into account the recoil of the molecules. This model leads 
to well-separated peaks at a given angle. The small angle elastic scattering, where the perfect 
gas model is no longer valid, should allow for the measurement of a coefficient of self diffusion. 
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I. Introduction. — De nombreux auteurs [1] ont 
donné un traitement complet de la collision d’un 
neutron avec une molécule d'hydrogène. Tenant 
compte de tous les degrés de liberté de celle-ci 
(niveaux de rotation, de vibration, états de spin) 
ils ont calculé dans le système du centre de gravité 
les amplitudes des diverses transitions. En som- 
mant sur l’angle de diffusion relatif, on pouvait 
aussi prévoir la section efficace totale associée à 
chaque transition. 

Si on envoie un faisceau de neutrons lents non 
polarisés sur de l'hydrogène liquide, c’est en réalité 
à un phénomène collectif que l’on aura affaire, et 
qui sera affecté par le recul ainsi que par les forces 
intermoléculaires. Grâce à l’effet du recul, la diffu- 
sion due même aux transitions élastiques sera très 
généralement inélastique. La question se pose de 
savoir si l’on pourra voir des niveaux de rotation 
dans le liquide. Il faut regarder en particulier si la 
diffusion pourra donner pour diverses transitions 
des raies résolues. Nous nous proposons de répondre 
à ces questions en calculant la section efficace 
différentielle, c’est-à-dire la répartition angulaire et 
énergétique des neutrons diffusés dans le système 
du laboratoire. Nous utiliserons le formalisme déve- 
loppé par R. J. Glauber et A. C. Zemach [2] et 
L. Van Hove [3] pour le traitement dans l’approxi- 
mation de Born de la diffusion par un système de 
particules. 


L’amplitude de diffusion par un proton se met 
sous la forme [2, 4] 


I1—2io 


_I1+1+2i0 
Fe 21 F1 


D SAS LCD 


a 


où /, i désignent le spin du proton, 
o le spin du neutron, 
a+ et a_ les amplitudes de diffusion correspondant 
aux noyaux composés { + 1/2 et I — 1/2. 
Nous appellerons par la suite Xe et 2 les carrés 
des amplitudes de diffusion « cohérente » et «inco- 
hérente ». 


CS PE se DA 
De — F RE CAE Dee a) — 0 07 410Æ2Darns 
TA) L 
Di — BI +1 (da )?—="1,586 barns. 


La coordonnée d’un proton sera mise sous la 


- Rire ; 
forme 7, + (—)*- où r, est la coordonnée du 


2 
centre de masse de la molécule d'indice /, A; la 
coordonnée relative des deux protons. &« est un 
indice égal à { ou 2, l'indice ! prend N valeurs. 
Nous négligerons dans À le degré de liberté 
correspondant aux mouvements de vibration de la 
molécule. C’est l’approximation du rotateur rigide, 
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approximation suffisante pour notre propos, car 
dans l’hydrogène liquide les molécules sont du 
point de vue de la vibration à l’état fondamental ; 
d’autre part l’énergie minimum d’excitation du 
premier état de vibration, de l’ordre de 0,5 eV, est 
très supérieure à celle des neutrons lents que nous 
considérons. E 

La coordonnée relative À représentera donc 
l’orientation de la molécule, à laquelle est associé 
le nombre quantique J. 

Il est raisonnable de faire les hypothèses sui- 
vantes sur le système diffuseur : 4) L’hamiltonien 
est indépendant des spins. 2) Il n’y a pas de corré- 
lation entre la position et l’orientation d’une molé- 
cule. 3) Il n’y a pas de corrélation entre orien- 
tations de molécules différentes. 


II. La fonction de corrélation. — Sous ces condi- 
tions la fonction de corrélation l'(7, t) caracté- 
ristique du système des V molécules se met sous 
la forme 


T = 4(2x)—$ (2N)—1 


N — —t — —— 

Do NdR exp - 2hr exp i Kri(0) exp 2 ÆKr(t) > 
1,j=1 

KRy(0) 


< al) exp — ie 


2 

DA 
p=1 

ag(t) exp i(—)8 0 > 


où se trouvent séparées les corrélations entre molé- 
cules différentes dues à leur mouvement de trans- 
lation, et les corrélations entre spins et orientations. 
La première moyenne 


fu = < exp — iKn(0) exp cKr;(t) > 
est caractéristique de la dynamique du liquide pro- 


prement dit et exige pour celle-ci un modèle. 
Quant à la deuxième, qui s’écrit 


aje(t) exp US > 


le modèle que nous avons choisi pour la molécule 
et les hypothèses (1) et (2) suffisent à la calculer. 
Ce terme g fera apparaître les effets d’interférence 
dus à la présence de deux protons dans la molécule 
ainsi que les possibilités de transition. Cette deu- 
xième moyenne doit porter sur les états internes 
des molécules caractérisés par les nombres quan- 
tiques associés aux À; et sur le spin du neutron 


pour tenir compte de la non polarisation du fais- 
ceau incident, 


INC 


La fonction de corrélation se met avec ces nota- 
tions sous la forme l 


(F4) — 4(2r)—$ (2N)—1 | 
fx exp — | 2 fu 85 + À fu |] = La + ls |R 
Hi 


où nous avons séparé les contributions provenant 


de molécules différentes (l,) et les contributions | 


provenant des protons d’une même molécule (1). 
l, joue le rôle d’une fonction d’autocorrélation | 
d’une molécule. | 


III. La diffusion « cohérente ». — Traitons | 
d’abord de la contribution du terme F4. 

Le terme g,,, pour L 7 s’écrit, compte tenu de | 
l’absence de corrélation d’orientation ou de spin 
entre deux molécules différentes, comme un produit 
de deux moyennes indépendantes, soit 


— 


. KR .KR 
DS no | & exp — ie + ds exp io ke 


os, est la projection sur Oz du spin du neutron. 
In > désigne un état d’une molécule, et px, son 
poids statistique. a, et a, sont les amplitudes de 
diffusion des deux protons d’une molécule. 

La quantité 


AT ir: 
A = QG exp — ie + 42 exp 


s'écrit 


Da RE 
PR SOS ES = 


An: ] cos 7 + 


2 re KR fes SO ER 
91 E 4 (a+ re Gien) | 53 cos À + LOU) sin ÀÀ | 


où S — 1, +1, est le spin de la molécule. 

L'état In, > se met sous la forme [JJ, SS, >, 
J'et J, décrivant la rotation de la molécule. 

Dans l’élément de matrice < n,|A|n, > la con- 
tribution de 6(1,—1,) sera nulle, cet opérateur 
n'ayant d’élément de matrice non nul qu'entre états 
de symétrie différente, c’est-à-dire de spin différent. 
La contribution à la moyenne totale de So sera 
aussi nulle, la trace des matrices de spin étant 
nulle, On obtient ainsi 


ue sin? (Ka/[2) 
C7 os er TELE 


où a = 0,75 À est la séparation des protons de la 
molécule d'hydrogène. On notera que g, ne dépend 


pas des poids statistiques des différents états pos- 
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sibles de la molécule d'hydrogène. La fonction T4 
s’écrit donc 


—- Sin? (Æa/]2) 
tKr (Kaf2)° 


DT). (2N) 14 >. fax exp 
> < exp — à Kry(0) exp & Krjlt) y. 
2] 


Le sens physique de T4 apparaît ainsi clairement. 
Représentant la corrélation entre protons de molé- 
cules différentes. nous y trouvons. 


a) Le facteur > . < EXp — iKri(0) exp iKr;(t) > 


qui correspond à He Aston entre positions + 
deux molécules. 


F - < sin? (Ka/2) 

b) Le facteur d’interférence Xe ECTDE qui 
fait intervenir la section efficace « cohérente » et 
qui vient simplement du fait que les protons sont 
liés par paires. 

Ces remarques suggèrent de définir la fonction 
de corrélation classique G4 associée aux mouve- 
ments des centres de masse des molécules. 


Galr, t) = (2x)—3 N—1 


f exp — ÎKr Z < exp — 5 KA(0) exp £K(4 >7. 


(1) 


La section efficace différentielle associée à Ty se 
met alors sous la forme 


sin? (Ka/2) 
Érva 12e (Ka/2)? 


f exp (Er ot) Gal, td) dr dt. (2) 


La section efficace associée à Ty est donc celle 
qu’on aurait pour un liquide monoatomique où 
chaque atome aurait une amplitude de diffusion 
égale à 4ücon Sin (Ka/2)/(Ka/f2) où Gcon — V2 est 
l'amplitude de diffusion « cohérente » du proton. 
La petite valeur de Gcon, alliée au facteur géomé- 
trique sin (Ka/2)/(Ka/2), rend la section efficace 
de diffusion « cohérente » très faible. 


IV. La contribution de I’, ou la diffusion « inco- 
hérente ». — Nous appellerons p, le poids statis- 
tique des molécules dans l’état J. En fait dans 
Phydrogène liquide n’existent pratiquement que les 
états J = 0 et J — 1 (para et ortho). L’énergie 
* associée au niveau J est 


SU TEA 
ss Eee — 0,015(1 /2) J(J 


+1)eV 


Le terme que nous avons appelé gy, qui repré- 


DIFFUSION DES NEUTRONS LENTS 


785 


sente les corrélations à l’intérieur de la molécule. 
se met aisément sous la forme 


eu = DPI ET ES 15 2 


e. 
2 IS J'JiS'SAI, 56, > [2 


4° g 


3 exp 7 Es - — Ey 


Rappelons que l’opérateur A se met sous la 
forme 


ae rt ne À KR 
: Tri s+ rat) es + 


2 EN D A RL ET OA 
ST (az — a_) [5 cos + io(i,— 1) sin 5 


Dans une transition conservant le spin de la 
molécule, seule la partie symétrique (soulignée) de 
l'opérateur donnera des éléments de matrice non 
nuls. 

Dans une transition changeant le spin, seule la 
partie antisymétrique interviendra. Dans une telle 
transition en particulier n’apparaîtra que ÿ’ampli- 
tude de diffusion «incohérente ». 

Nous appellerons gy->y la quantité g4 corres- 
pondant à une transition J + J'. Le calcul montre 
qu’elle s’écrit suivant le cas, pour une transition 
conservant le spin : 


PJ URSS 


Eee “Xe + 


BIT! = + 4/88(S +1)2i] 


ñ 


on JT Loos ÆR jrs, HE 


Pour une transition changeant le spin : 


ï .Eyr — Ey SUIS. = 
DT  - exp > : fi ] & Di 
AY Le D KR 2 
PA < J'Ji|sin | > 
Posons 
Eyr 


mes IGN, : 
JT = EXP E——— tjr (A). 


Définissons d’autre part la fonction d’auto-corré- 
lation des centres de masse, qui est une caracté- 
ristique de la dynamique du liquide, 


Gr, 1) = (2x) 


Jarexp — i kr < exp — à Kry(o) exp iKrilt) 7. (3) 
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La section eflicace différentielle associée à la tran- 
sition J — J'se met alors aisément sous la forme 


ke 
[— 14 T D S LA 
jo de — *NUr) 19 (K) 


fepifr- (eo) Jar aa de (0 


ñ 


Ce résultat appelle deux remarques : 

a) L'apparition dans la transformée de Fourier 
de © —(Ey — Ej)/,; comme argument du temps, 
au lieu de « provient de ce qu’en plus des effets 
de la dynamique du liquide proprement dit inter- 
vient aussi dans le bilan énergétique l’énergie asso- 
ciée à la transition. 

De façon plus précise appelons V, la vitesse ini- 
tiale de la molécule. Compte tenu de la conser- 
vation de l’impulsion, l’équation de conservation 
de l’énergie s'écrit 


Onre n K)? 


ho — Eyr E | 
"| Va &m 


2 
mvVè 


ou 


ho (Es Ey) = LH 


Le 2° nombre étant en moyenne nul, on peut 
écrire une équation moyenne de conservation de 
Ci : 
l’énergie 
72 K°? 

SE tu (5) 
km 


re) (Ey, Ey) 
qui détermine l’emplacement en énergie du pic 
diffusé. Cette équation tient compte comme on le 
voit de l’énergie de recul de la molécule qui ne 
peut jamais être négligée dans le cas de la molécule 
d'hydrogène. Le recul peut en effet être négligé 
(ouaV) sen Te LT où mret Te mMVetoT 
désignent respectivement la masse et la tempé- 
rature du neutron et de l’atome diffuseur. Ce rap- 
port, égal 1c1 à T,/2T (T = 20 0K)est toujours > 1 
pour des neutrons de température > 20 0K. 

b) La quantité /5_,y(K) qui joue le rôle du carré 
d’une amplitude de diffusion (dépendante de 
l’angle) est déterminée par la transition considérée. 


081 Bamns/s'éradian. 


%4 


— 


PR nn re 
90 Angle en degrés 180 


F16. 1. — do/d{2 par molécule pour la transition 0 — 0 
pour diverses longueurs d’onde initiales données en 
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barns/stéracian. 


90 Angle en obgrés 180 


| 
F1G. 2. — do/dQ par molécule pour la transition 1 — 1 
pour diverses longueurs d’onde initiales données en À. 


Darns /Stéradian . 


08 


à 


À = 1141 À 


90 Angle en degres 180 
F1G. 3. — do /dQ par molécule pour la transition 4 — 0. 
barns/shéradian. 
42 
8 
4 
0 —————— ;— 
30 Angle en degres ‘80 
F16. 4. — do fdQ par molécule pour la transition 0 — 1. 
À = 1,141 


Son ordre de grandeur varie fortement d’une tran- 
sition à l’autre. 


Voici son expression pour quelques transitions, 
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16 
barns/stéradian. 


0B 


a ———— 


F1c. 5. — do /dQ par molécule pour la transition 1 — 2. 
ko = 1,141 À. 


les seules possibles avec des neutrons de longueur 
d'onde > { À 


Îo—0 — 4Po 2e 16 
foi = 4PilDe + 2/8 Di] [58 + 272] 
fo = Pi Z:/3 oh 
foi = 4po 32 j? 
fi=>2 =" 4p DE (275 + 85). 
Les 7 désignent les fonctions de Bessel sphériques 


prises pour la valeur Xa/2. 


V. Expression de la section efficace dans un 
modèle de gaz parfait. — Revenons à la formule (4) 
qui donne la section efficace différentielle associée 
à une transition 


d?6. 
dQ de 


k 
NOT) FT Le ES 


fexpif#r (o Er : 2) ( Gr, t) dr de, 


Si on néglige l’effet des forces intermoléculaires, 
on peut prendre pour G;(r, t) la fonction d’auto- 
corrélation du gaz parfait 

27 (kg Ti— 1%)? 
2m | ae 21 


2mr? 


(kp Tt— 1à) 


GF,0 = | 6) 
(CL TIT): 

L’approximation ainsi faite sera toujours très 
bonne aux temps courts, c’est-à-dire pour 


/ 


E y ii Ey 


Rte grand. L’équation (5) montre que 


grâce au recul il en est toujours ainsi sauf dans le 
cas des transitions élastiques aux petits angles. 
Nous reviendrons sur ce cas particulier, On ne 
peut dans lexpression (6) négliger # sous peine de 
négliger le recul (cf. V). On obtient avec ce modèle 
la formule de section efficace suivante par molécule 


d?6 ï i mn 2 k fy>gJ'\K) 
dQ de Tr kg T #ko K 


Uro) (E y 
EXP — m — 


$ 1H 
3) En | 


KT EE 2 


Où l’on voit apparaître comme prévu la conser- 
vation de l’énergie. 


F1G. 6. — Surface de diffusion associées aux diverses tran- 
sitions pour une longueur d’onde initiale 14 — 1,141 À. 
Valeurs de la longueur { 1 —> 0 Ny = A1,60À 
d’onde finale À, 0—0, 1—>1 À, — 2,12 À 
pour l’angle 0 = xf2 (0 —1 = 4 000 


VI. Les surfaces de diffusion associées aux tran- 
sitions. — [’équation (5) de conservation de l’éner- 
gie, qui donne l’emplacement de la raie diffusée, 
s'écrit 

 K2— k? 


n2 = 
2m 


DRE 


km 


be Lo) 


où À — k-— ko, ko et k vecteurs d'onde du neu- 
tron incident et diffusé. Dans l’espace des k, le 
lieu des extrémités des vecteurs # est donc une 
sphère. Le centre de ces sphères est le point C 
(voir figure) défini par OC = k,/3 (O origine des 
vecteurs). Le rayon de la sphère correspondant à la 
transition J — J'est 


Fe: Ê : 


Sp E] 1 
E, 
en appelant £, l'énergie du neutron incident. 
La sphère À de rayon 2/3 k, correspond aux 
transitions élastiques ; Z, et X, représentent des 
transitions où la molécule est excitée à un niveau 
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supérieur. », représente une transition où la molé- 
cule se désexcite. Dans les deux premiers cas le 
neutron perd toujours de l’énergie. Dans le troi- 
sième cas il n’en gegne pas toujours. 


Pro 


A un angle donné les raies correspondant aux 
diverses transitions seront nettement séparées en 
énergie si leur largeur est faible vis-à-vis de leur 
espacement. Il faut faire exception pour les raies 


élastiques 1 — 1 et 00 qui seront toujours 


confondues. 


VII. Largeur des raies. — Avec les hypothèses 
que nous avons faites, la largeur de la raie diffusée 
provient uniquement de l’agitation thermique des 
molécules. 

Dans l’expression (7) de la section efficace, la 
largeur est essentiellement imputable à l’exponen- 
tielle, qui s’écrit 


K? 2 


où À est le rayon de la sphère de diffusion, et £, 
l’énergie du neutron incident. 

Dans tous les cas la largeur est minimum aux 
petits angles. On se rend compte qu’elle a pour 
expression quand elle est faible 


Kk, 
3k— k, 3 


TT 1/2 
dk (3 =) 


E 


Nesit 


La largeur en énergie correspondante est 


kp T 


1/2 
E—uar) Î(æ— 1) (1 + 2a)| 


9 , 


où nous avons posé 


3/2 5e 


Il est facile de voir que cette largeur en énergie 
décroit dans tous les cas quand l’énergie inci- 
dente E, croît. Ainsi pour Æ, — 720 ©K on trouve 
respectivement pour les transitions 0 — 1 et 1 0 
des largeurs de l’ordre de 20 et 13 ©K. Pour une 
énergie incidente de 20 °K et la transition 1 — 0 
on trouve une largeur de 30 0K, l’énergie finale 
étant alors de l’ordre de 180 °K. 

La largeur pour les transitions inélastiques varie 
donc peu autour de 20 ©K. Le gain ou la perte 
d’énergie dû aux transitions étant beaucoup plus 
grand, les raies sont en général bien séparées. Leur 
largeur est cependant trop forte pour pouvoir 
espérer voir la largeur propre des niveaux de rota- 
tion dans le liquide. 


Cas particulier des transitions élastiques aux 
petits angles. — Comme nous l’avons remarqué 
plus haut, le modèle de gaz parfait est incorrect 
pour expliquer la diffusion élastique aux petits 
angles. Si l’angle est suffisamment petit, l’échelle 
des temps est grande et la fonction d’autocorré- 
lation peut être déduite d’un modèle de self diffu- 
sion. La condition à laquelle on est conduit pour 
l'angle est en fait Xd x (cf. V), qui s’écrit encore 


10 > 0d 


où À, est la longueur d’onde initiale, 0 l’angle et d 
une distance intermoléculaire moyenne. Avec des 
neutrons de l’ordre de 5 À et des angles inférieurs 
à 100 il doit donc être facile de mesurer le coef- 
ficient de self diffusion de l’hydrogène. 
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SOUS-STRUCTURES DANS LES CRISTAUX DE BROMURE D'ARGENT 


Par Henri CHATEAU et ANNE DE CUGNAC, 
Laboratoires de Recherches de la Société Kodak Pathé (1) 


Résumé. — On a étudié par diffraction de rayons X monochromatiques des monocristaux de 


bromure d’argent de forme très régulière, 


D’après les diagrammes, chaque cristal est constitué par un réseau régulier au sein duquel 
existent de très nombreux petits cristallites orientés par rapport aux axes [110] du cristal de base. 
Les angles de rotation des cristallites autour de ces axes sont généralement inférieurs à 10°, mais 


peuvent atteindre 260. 


| Les sous-structures se forment à l’obscurité pendant la croissance des cristaux ou ercore, sous 
l'action de la lumière pendant le noircissement direct. 


Abstract. — Several monocrystals of silver bromide were analyzed by monochromatic X-ray 


diffraction. 


The shape of these crystals is very regular, nevertheless the diffraction patterns 


show the presence of a great number of substructures in every crystal. 

According to the patterns, each crystal is constituted by a regular lattice in which there are 
small crystallites. These crystallites are oriented around the 110 axis of the perfect crystal and 
the angles of rotation are generally less than 10°, but in some cases, they reach 26°. 

The substructures appear in darkness during the growth of crystals or, also, whilst the light 


causes blackening of crystals. 


La diffraction de rayons X effectuée sur de nom- 
breux monocristaux de BrAg préparés par voie 
humide a donné des diagrammes de Laue d’aspect 
inattendu. À côté des taches normales présentant 
un astérisme considérable, on observe des taches 
propres à une structure fibreuse. Celles-ci révélant 
lPexistence probable de sous-structures ou cristal- 
lites, nous avons approfondi l’étude de ces cristaux. 

Les cristaux ont été obtenus à partir de diverses 
solutions de BrAg dans des solvants appropriés : 

1° BrAg est dissous dans l’ammoniaque con- 
centré bouillant. Lorsque la solution est saturée on 
transvase dans un bécher maintenu à 20 0C. Il se 
forme en quelques minutes des cristaux de BrAg. 
Ceux formés à la surface du liquide sont plus gros 
et plus homogènes en taille que ceux formés au 
sein de la solution. 

20 De l’acide bromhydrique est saturé en BrAg 
à 70 0C. Par refroidissement de gros cristaux se 
déposent. 

30 On porte à 70 0C de la diméthylformamide 
préalablement saturée à 25 © en BrK. On ajoute un 
excès de BrAg. Après filtration, la solution refroi- 
die à 250 donne de beaux cristaux. 

Ces trois techniques permettent l’obtention de 
cristaux tabulaires dont la taille est généralement 
comprise entre 100 et 700 w. L’épaisseur peut 
atteindre 20 u, mais elle est souvent de l’ordre 
de 5 pu. 

Le contour des faces à grand développement est 
presque toujours très simple : hexagone régulier ou 
semi-régulier, triangle équilatéral. 


(*) Communication n° 2079 V. 


Technique d’analyse aux rayons X. — Le cristal 
est fixé sur un support mince d’acétate de cellulose 
par une bande collante ou par une goutte d’une 
solution de gélatine ou de graisse de silicone, Il est 
ensuite orienté de telle façon que les faces de grand 
développement soient perpendiculaires au faisceau 
de rayons X, la position des axes cristallins étant 
repérée. 

Le générateur est un « Dictyobloc » à tube scellé. 
L’anticathode est en cuivre et le rayonnement émis, 
bien que filtré par du nickel, possède un résidu 
polychromatique. La tension d’alimentation est de 
45 kV et le débit de 6 mA. 

La divergence des rayons X est réduite à 10 par 
une double collimation. 

Dans ces conditions, le temps d'exposition varie 
de 2 à 12 heures pour obtenir les taches de diffrac- 
tion de texture. 


Repérage des axes. — Pour repérer les axes 
cristallographiques, on a appliqué la technique du 
cristal tournant. Dans tous les cas, on a vérifié que 
les indices sont bien ceux supposés en admettant 
que les faces de grand développement du cristal 
correspondent à des plans { 111 #. 

La photographie reproduite figure 1, correspond 
au diagramme d’un cristal de BrAg de 300 w 
tournant autour d’un axe [110]. On remarquera 
que les taches de diffraction sont allongées suivant 
l’arc Debye-Scherrer passant par ce point. Nous 
verrons l'importance de cette observation. 


Aspect des diagrammes de diffraction. — 
L'exemple type représenté figure 2 a été obtenu 
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Hypothèses d’axes de fibre [110]. Caleul du 
nombre et de la position des taches de diffraction. — |}: 
Chaque désorientation peut être considérée comme |} 
une partie du cristal initial parfait ayant subi une |h 
rotation autour d’un axe (analogue à un axe de |hn 
fibre). On peut rechercher l’indice des divers axes |hh 


avec un cristal de 250 un (fig. 3). Les taches de 
diffraction présentent une symétrie ternaire ou 
sénaire. Elles sont situées sur des cercles concen- 
triques qui correspondent aux anneaux Debye- 
Scherrer d’une poudre de BrAg. 
Des taches présentant la même symétrie appa- 
raissent quelquefois entre les anneaux 200 et 220 


du BrAg. Ces taches sont dues à de l'argent pro- 


venant d’une photolyse et correspondent à la dif- 
fraction sur les plans { 111 }. 

Lorsque l’axe vertical du cristal est un axe [2111], 
le diagramme a l’aspect de la figure 2. L’axe ver- 
tical du diagramme traverse une tache sur chacun 
des anneaux 400, 222 ou 200. 

La grandeur et la forme des taches de diffraction 
varient légèrement d’un cristal à l’autre. Les taches 
sont parfois ponctuelles mais généralement elles 
forment un arc sur l’anneau Debye-Scherrer corres- 
pondant, l’angle de cet arc étant toujours inférieur 
à 100 et le plus souvent de l’ordre de 2 à 4°. 
Lorsqu'une tache constitue un are de 4 à 100, on 
peut souvent la décomposer en plusieurs points 
d’intensités relatives différentes (fig. 4), ce qui est 
l’indice de la présence de sous-structures dans le 
cristal. 

Les plus gros cristaux (> 400 u) donnent sou- 
vent des taches de Laue avec astérisme (fig. 5). 
Leur allongement considérable est propre à des 
cristaux imparfaits. 


Interprétation des diagrammes de rayons X. — 
Pour un cristal parfait en position fixe dans un 
faisceau de rayons X monochromatiques, l’orien- 
tation satisfaisant à la relation de Bragg n’est 
réalisée qu’accidentellement. 

Les cristaux analysés ayant toujours un axe [111] 
parallèle aux rayons X, le calcul montre qu'il est 
impossible avec un cristal parfait d'obtenir des 
taches de diffraction sauf, peut-être pour le plan 
(400) (la divergence des rayons X étant de 10 et 
l'orientation des cristaux tabulaires ne pouvant 
être réalisée à mieux de 1 à 29 près). 

Pour interpréter les diagrammes de diffraction, 
il faut admettre que le cristal présente des zones de 
désorientation pour lesquelles les plans cristallins 
ont subi des déviations angulaires. 

Nous avons vérifié que les rayons X sont inca- 
pables de produire le type de désorganisation 
observé, même après plusieurs dizaines d'heures 
d'exposition sur le diffractomètre. Cette obser- 
vation est en accord avec les résultats de 
R. H. Herz [1] et de’ C. R. Berry [2-3]. 

Le petit nombre de taches observées et leur 
symétrie indiquent que les désorientations ne sont 
pas quelconques. 

Les diagrammes obtenus sont semblables aux 
diagrammes de textures métalliques. Ceci suggère 
qu’un axe donné et commun à toutes les zones 
désorientées garde sa direction initiale. 


IN] 


de rotation possibles en considérant les angles que 


font les droites passant par les taches et le centre (l 
du diagramme avec l’axe vertical du film. Cette 1 


méthode d’analyse est classique dans l’interpré- 


tation d’un diagramme de texture. Mais si, d’une | 


manière générale, il n’y a qu’un ou deux axes 


de fibre, dans le cas des cristaux de BrAg on peut | 


prévoir plusieurs «axes de fibre ». 


Soit à l’angle que fait la droite joignant une | 


tache au centre du diagramme et la verticale du 
film. M. Polanyi a montré (4) que 
cos à — cos p/cos 0 

où p est l’angle que. fait le plan cristallin de la 
tache considérée avec l’axe de fibre et 0 l’angle de 
Bragg. 

Si le cristal est orienté de telle façon qu’un 
axe [211] soit vertical, on obtient le diagramme 
schématisé figure 6. Une tache de chacun des 


Axe 211 Vertical 


400 
ee 
200 220 
eee / 200 \ gee 
im, sm 
0 280 
200 | | 
200 \ / 200 
k 
in 
220 7 280 
400 \ / + 
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Fic. 6. — Schéma de diagramme plan de diffraction de 
rayons X d’un monocristal de BrAg. 


plans (400), (200) et (222) est située sur la verti- 
cale. Il n’y a que deux valeurs expérimentales 
de à : 0 et 600 pour les plans cités, 30 et 90° pour 
le plan (220). 

Le fait d'obtenir à — 0 pour des plans d'indices 
différents indique qu’il y a plus d’un axe de fibre 
et probablement trois en raison de la symétrie de 
toutes les taches. | 


F16. 1. — Cliché de diffraction d’un cristal de BrAg 
tournant autour d’un axe [110]. 


Fic. 3. — Monocristal de BrAg. 


FiG. 2. — Diagramme de diffraction 
d’un monocristal de BrAg. 


F1G. 4. — Agrandissement d’une tache 
de diffraction de BrAg. 


Fic. 5. — Cliché de diffraction 
montrant des taches de Laue avec astérisme. 
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F1G. 7, — Diagramme de poudre du BrAg. F16. 8. — Cristal ayant subi une corrosion prolongée. 


Fi, 9. — Diffraction d’un cristal peu déformé F16. 10. — Diffraction du cristal de la figure 9 
avant exposition à la lumière. exposé 4 heures à la lumière solaire. 


* No! 


D'autre part, les angles à ont souvent les valeurs 
30 ou 600, à mieux de 10 près. Or, il existe trois 
axes [110] dans un plan (111) faisant deux à deux 
un angle de 600. Il en est de même pour les 
axes [2111]. En outre, l’angle entre un axe [110] °t 
un axe [211] est égal à 300 ou 900. 

Ces remarques montrent la’ possibilité d’avoir 
trois axes de fibre ayant les mêmes indices : [110] 
ou [211], ces axes étant situés dans un même 
plan (111). 

Avec des axes de fibres [211], chaque tache 200 
conduit à une valeur p de 66° qui, reportée dans la 
formule de Polanyi, conduit à 8 — 650, Comme 
l'angle expérimental est presque toujours égal à 600, 
il y à incompatibilité. Cette conclusion est con- 
firmée par la construction d'Ewald qui permet le 
calcul exact de à. 

Au contraire, si l’on admet que les axes de fibre 
sont les trois axes [110] situés dans le plan du 
cristal, on prévoit très exactement l’emplacement 
sur le cliché de toutes les taches observées. 

Mais la théorie prévoit beaucoup plus de taches 
qu’il n’en est observé. Le problème est analogue 
lorsqu'on considère une tôle laminée, certaines 
taches manquant dans les diagrammes. Par ana- 
logie, on peut supposer que dans le BrAg les 
plans Ÿ 111 } parallèles aux faces de grand dévelop- 
pement peuvent s’écarter de leur position normale 
pour créer une désorientation, mais très diffci- 
lement au delà d’un angle limite. 

Connaissant les axes de fibres, il est possible de 
prévoir l'emplacement et le nombre des taches de 
diffraction en fonction de l’angle &« de rotation. 

Pour faciliter l'exposé, on suppose un axe de 
fibre [110] vertical. L’axe [111] perpendiculaire au 
cristal est confondu avec les rayons X. Le problème 
de la position théorique des taches est alors résolu 
graphiquement à l’aide de la construction d’Ewald. 
Pour les tracés, on utilise les données suivantes : 


constante du réseau : 5,775 À 
longueur d’onde Æo : 1,541 À (2/3 Ka; + 1/3 Ka). 


Soit y l’angle de la projection d’un rayon dif- 
fracté sur le plan horizontal avec le rayon incident. 
On a déterminé y à mieux de 1/20 près. 

Connaissant y et 0 on peut calculer à par la 
relation : tg y — tg 20 sin à. Il est possible de 
contrôler rapidement cette façon de procéder en 
caleulant ?, angle entre l’axe de fibre et la normale 
au plan étudié, par la relation : cos 9 — cos 0 cos à. 

La valeur obtenue doit coïncider avec une des 
valeurs possibles déterminés par la théorie. 

Les résultats obtenus ont permis d'établir le 
tableau I qui considère uniquement les trois axes 
de fibres [110] situés dans le plan du cristal. 

Les angles à qui conduisent aux taches de dif- 
fraction observées sont en italique. On voit qu'à 
tous ces angles, il est possible d’adjoindre un angle 
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de rotation & toujours inférieur ou égal à 260 (ces 
valeurs sont en italique). Lorsque « est plus grand, 
avec l’appareillage utilisé, 11 n’y à jamais de taches 
de diffraction visible, 


TABLEAU I 
PLAN Y p 0 “ 
THÉOR. EXP. ? ANEPE 
DE ROTATION 

200 022204: 60, CAGED US S80 780 
319 900 909 90° 19° 519 
220 3991 07 60° 60035’ 97019 260 459 820 
1 559 900 900 909 320 770 
999 299300085016 2M55095 230 189 880 
a 550 90° 900 900 8° 89 62° 
400 50° 459 LE CHA 349 70 760 
6590 90° 900 90° 39 679 
111 459 350106051085 05304% 07110 599 
26030” 90° 83° 90° 69° 349 770 
2208051050 8510207 48° SAS 379 770 
311 48° 65050’ 65008’ 61040 6° 459 659 770 
52930 900 89° 90° 479 399 86° 


En résumé : on prévoit les taches de difjraction 
observées en admettant que le cristal possède des zones 
désorientés d’un angle inférieur à 30° autour d'un 
axe [110]. Le rôle des trois directions [110] non 
situées dans le plan du cristal n'apparaît pas dans 
les diagrammes. Ces sous-structures peuvent être 
constituées soit par des blocs mosaïque, soit par des 
portions de cristal déformées. Dans ce dernier cas, 
il y a continuité entre Le cristal de base et les régions 
désorientées. 

Les diagrammes de rayons X ne permettent pas 
de choisir entre ces deux représentations. Pour 
simplifier cet exposé, on utilisera au cours de linter- 
prétation des résultats le terme « cristallites » pour 
désigner les régions désorienteës sans préjuger de leur 
nature. 

Les cristallites doivent être nombreux car les 
taches de diffraction sont relativement intenses. En 
outre, ils sont certainement assez petits pour que Île 
cristal garde sa cohésion et sa forme régulière 
macroscopique. 


Interprétation de l’intensité relative des taches 
de diffraction. -— Comme toutes les désorientations 
ne participent pas à la réflexion sélective de Bragg 
et que la probabilité d’un angle de rotation éleve 
doit être plus petite que celle d’un angle faible, les 
taches pour lesquelles l'angle &« est petit doivent 
avoir une intensité relativeinent plus grande. Mais 
il faut tenir compte de la géomécrie de la diffrac- 
tion et du facteur de structure des plans eristallins. 

La figure 7 représente un diagramme Debye- 
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Scherrer d’une poudre de BrAg. Les intensités des 
anneaux décroissent quand le diamètre grandit. 
Au contraire, avec les monocristaux les taches 
400 sont intenses alors que les taches 220, par 
exemple, sont très faibles. Mais dans le tableau I, 
l’angle « est de 30 pour les taches 400 alors qu’il 


est de 26° pour les taches 220. Cette inversion des 


intensités s'explique donc aisément. 

D'une façon générale, pour distinguer une tache 
sur le cliché, il est nécessaire d’avoir un angle « 
d'autant plus petit que l’anneau Debye-Scherrer 
donné par la poudre correspondante a une intensité 
plus faible. 

Pour une durée d’exposition aux rayons X ne 
dépassant pas cinq heures, beaucoup de cristaux 
donnent des taches bien visibles. Ce fait est remar- 
quable car seuls interviennent dans la diffraction 
les cristallites désorientés. Cependant, il faut noter 
que chaque tache de diffraction est formée par les 
rayons diffractés provenant de tous les cristallites 
pour lesquels la rotation a amené le plan considéré 
sous l’angle de Bragg. Il y a convergence au même 
point, au contraire d’un diagramme de poudre où 
chaque grain donne une tache distincte. Chaque 
tache de diffraction étant due à la présence de 
plusieurs cristallites, il n’est pas possible d'évaluer 
le nombre de cristallites par cristal. 


Relation entre les sous-structures des grains de 
BrAg et la dimension des taches des diagrammes 
de Laue ou des diagrammes de eristal tournant. — 
Le cristal contenant des cristallites désorientés, on 
doit observer un diagramme de Laue avec asté- 
risme et la longueur de l’astérisme doit être direc- 
tement fonction de la dispersion angulaire de déso- 
rientation. Les extrémités de la tache d’astérisme 
sont faibles car elles correspondent aux cristallites 
les plus désorientés donc en nombre moindre. 

Pour plusieurs clichés, on a trouvé des longueurs 
d’astérisme atteignant H 50, Ce résultat est infé- 
rieur à la limite déterminée précédemment, mais il 
correspond aux cas les plus fréquents de désorien- 
tation. 

Les diagrammes de cristal tournant permettent 
une bien meilleure appréciation de la désorien- 
tation des cristallites. Dans le cas où les taches de 
diffraction sont sur la ligne équateur, on montre 
qu’une désorientation « autour d’un axe [110] per- 
pendiculaire à la verticale provoque un dépla- 
cement angulaire 8 vertical de la tache par rapport 
au centre de la chambre tel que l’on ait : 


Nr ERREURS 
sin É) HS () .sin 0 


où 0 est l’angle de Bragg. 

On calcule que «& est compris entre 40 et 159 mais, 
pour des cristaux très perturbés, on a noté des 
valeurs atteignant — 230. On obtient dans ce cas un 


N°11 |b 


recoupement excellent avec les diagrammes de 
cristaux fixes. 

Les taches de Laue avec astérisme sont surtout 
visibles avec les plus gros cristaux (500 à 700 u), 


parce que le volume de matière irradiée est suffi- JM 


samment grand. Il est possible d’autre part que les 


gros cristaux tabulaires soient, plus facilement que | 


les petits, sujets à des désorientations car ils sont 
plus flexibles. 


Or, l’existence de déformations plastiques dans | 
un cristal de bromure d’argent a été récemment |} 


observée par R. H. Herz [5] qui a montré que | 
l’angle de désorientation pouvait atteindre 100. | 


Influence de la présence d’un dépôt microeris- | 
tallin de BrAg à la surface des monocristaux de |h 


BrAg. — On pourrait expliquer l’apparition des | 
taches Debye observées, en supposant que lors de |k 


la préparation des cristaux, il se forme un dépôt | 
orienté de microcristaux sur leurs faces à grand | 
développement. Cependant, ce dépôt est peu pro- | 
bable car il faudrait supposer que les microcristaux | 
ne soient pas à plat sur le gros cristal support et 
par suite qu’ils ne soient pas en épitaxie simple. 

Toutefois, pour réfuter cette hypothèse, nous 
avons effectué deux séries d’expériences. 

x) On a corrodé les cristaux par une solution 
diluée d’hyposulfite de sodium ou de bromure de 
potassium. Le lavage a été effectué de façon à 
éviter toute précipitation de BrAg par destruction 
du complexe argentique formé. Le cristal a été 
attaqué plusieurs fois jusqu’à corrosion apparente 
du cristal (fig. 8). Dans ces conditions, on peut 
admettre que la surface du cristal a complètement 
été débarrassée des microcristaux superficiels. Or, 
les diagrammes de diffraction ne révèlent prati- 
quement aucune différence. Il est prouvé que des 
cristallites existent dans le cristal et qu’ils ne sont 
pas principalement localisés à la surface. 

B) On a déposé des microcristaux sur un mono- 
cristal ne donnant pas de taches d'orientation | 
Debye-Scherrer. On a recouvert le cristal d’une 
solution ammoniacale saturée en BrAg, puis ajouté | 
un peu d’eau pour provoquer la précipitation de 
microscristaux de BrAg. Le diagramme de dif- 
fraction des cristaux ainsi traités montre qu'il ya | 
bien eu dépôt de petits cristaux mais sans orien- 
tation privilégiée. 

Ces expériences prouvent qu’il n’y a pas inter- 
férence dans la diffraction entre les désorgani- | 
sations orientées et un hypothétique dépôt micro- | 
cristallin orienté sur les gros cristaux. | 


Désorganisations orientées provoquées par pho- 
tolyse. — Les clichés de diffraction témoignent de 
la présence indubitable de cristallites dans des cris- 
taux préparés à l’obscurité. On a en outre observé 
que les dislocations analogues peuvent être pro- 
duites par la lumière (mais non par les rayons X) 


. 
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Des cristaux sans sous-structures apparentes 


- (fig. 9) ont été exposés à la lumière solaire et des 


chchés de diffraction furent pris après différents 
temps d'exposition. Au bout d’une heure, au dia- 


… gramme de Laue, se superposent de très faibles 


taches Debye-Scherrer d’orientation [110]. Puis 
ces taches s’intensifient pour des durées d’expo- 
sition plus longues. Au bout de six heures d’expo- 
sition, on obtient un intense diagramme de dislo- 


- cations orientées suivant [110], tandis que les 


taches de Laue ont presque complètement disparu 
(fig. 10). Sur ces diagrammes, on note aussi l’appa- 
rition des anneaux Debye-Scherrer de l'argent. 


Diseussion et conclusions. — Nous avons établi 
que la principale désorientation des cristallites 


— existant dans un monocristal de BrAg est une rota- 


tion autour d’un axe dont l'indice est toujours 110. 
Ce type de désorientation peut provenir du mode 
de croissance des cristaux de BrAg. Il est, en effet, 
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intéressant de remarquer que les arêtes d’un cristal 
de BrAg sont constituées par des rangées [110]. 

Lorsque le cristal grossit, l’arête initialement 
siège d’une désorientation cristalline, se trouve à 
l’intérieur du cristal. Ce schéma, qui n’est qu’une 
hypothèse, peut cependant convenablement expli- 
quer la présence des sous-structures observées à 
l’intérieur des cristaux de BrAg. 

Cependant, la formation du même type de sous- 
structure à été observée sur des cristaux exposés 
à la lumière solaire. [1 est probable que les ran- 
gées [110] du cristal de bromure d’argent sont rela- 
tivement plus résistantes que les autres et que les 
zones de rupture sont toujours les mêmes quelle 
qu’en soit la cause. 

L'existence des sous-structures orientées dans le 
bromure d’argent est en accord avec les hypothèses 
de quelques auteurs et notamment de C. R. Berry 
et R. L. Griffith sur le mécanisme du noircis- 
sement direct à la lumière. 


Manuscrit reçu le 5 juillet 1960. 
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MESURE DES INDICES DE RÉFRACTION DANS L'INFRAROUGE 


Par Mme JoserrEe VINCENT-GEISSE et M. Jean LECOMTE, 
Laboratoire de Recherches Physiques, Sorbonne. 


Résumé. — Les mesures d'indices de réfraction dans l’infrarouge ont été très nombreuses ces 
dernières années. Après avoir rappelé l’intérêt qu’elles présentent dans différents domaines, nous 
décrivons les méthodes utilisées, en insistant sur les plus nouvelles, et nous donnons une biblio- 
graphie des travaux effectués récemment sur cette question. 


Abstract. — There have been numerous measurements of refractive indices in the infrared in 


recent years. 


After recalling the interest which they present in different domains, we have des- 


cribed the methods utilized, principally the newer ones, and compiled a bibliography of the work 


recently carried out on this subject. 


Les mesures d'indices dans l’infrarouge, bien que 
très anciennes, étaient relativement peu nombreuses 
jusqu’à ces dernières années, et se présentaient souvent 
comme une simple extension des mesures effectuées 
dans la région visible du spectre. 

L'un d’entre nous avait publié en 1955 [42] une 
mise au point des méthodes utilisées et des résultats 
obtenus et concluait en ces termes : « La réfracto- 
métrie dans l’infrarouge reste encore un vaste champ 
ouvert aux déterminations expérimentales et aux véri- 
fications des théories. » 

Les mesures d'indices connaissent, depuis cette 
époque, précisément, un regain de faveur, les publi- 
cations relatives à ce sujet, datant seulement des cinq 
dernières années, étant aussi nombreuses que celles qui 
avaient paru auparavant. 

Nous avons donc pensé qu’il serait utile de reprendre 
la question, d'examiner les raisons qui motivent le 
renouveau de son intérêt et de décrire les méthodes 
nouvellement utilisées comme les résultats obtenus. 


Introduction. Intérêt des mesures dans l’infrarouge 
d’indices de réfraction. — 1° Encore récemment les 
mesures d'indices dans l’infrarouge présentaient un 
intérêt principalement à deux points de vue : 

a) Recherche dans un domaine étendu de longueurs 
d’onde, de cristaux transparents et présentant la plus 
grande dispersion possible, C’est dans ce but que l’on a 
sélectionné successivement, par exemple, le bromo- 
iodure de thallium, le bromure, puis liodure de césium. 

b) Vérification des différentes formules de dispersion 
proposées, comme celles de Cauchy ou de Ketteler- 
Helmholtz, et prévision des longueurs d’onde propres 
d'absorption. 

29 Ces dernières années, différents travaux relatifs 
à d’autres sujets, ont conduit à exécuter de nouvelles 
mesures d'indices de réfraction dans l’infrarouge et 


l'intérêt porté à ces mesures s’est considérablement 
élargi et diversifié. 

a) La détermination de l'intensité des bandes 
d'absorption de corps en solution exige la connaissance 
de l’indice du solvant : on a été ainsi amené à mesurer 
l'indice des solvants usuels dans l’infrarouge. 

b) La théorie des vibrations des ions, à l’intérieur 
d’un cristal ionique, fait intervenir des fréquences 
propres de vibration ; celles-ci sont difficiles à obtenir 
par les méthodes usuelles de transmission, à cause de 
l’absorption considérable à l’intérieur de ces bandes. 
La détermination de l’indice de réfraction, jointe à celle 
de l’indice d’absorption, fournit, avec précision, les 
fréquences cherchées. Quelques mesures d’indices ont 
été effectuées dans ce but. 

c) D’une manière plus générale, intensité de bande 
et indice de réfraction au voisinage d’une bande sont 
liés dans les théories de la dispersion. Les mesures 
d'indices de réfraction sont utiles, d’une part, pour 
vérifier ces théories, et, d’autre part, pour remplacer 
les mesures d’intensités partout où ces dernières se 
montrent difficiles, c’est-à-dire principalement dans le 
cas des bandes très fortes. 

d) Si nous passons aux métaux, la théorie de Drude 
relie les constantes optiques à un certain nombre de 
grandeurs, conductibilité, concentration des électrons 
libres, temps de relaxation ; leur mesure directe est. 
donc également intéressante. 

e) Enfin, on a assisté, ces dernières années, à un 
développement considérable de l’utilisation des semi- 
conducteurs. Là encore la mesure des constantes 
optiques présente un intérêt théorique et pratique de 
premier ordre. 

f) Nous citerons encore pour mémoire, la construc- 
tion de filtres interférentiels ou de couches antiréflé- 
chissantes, qui exige la connaissance des constantes 
optiques d’un grand nombre de corps. 
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mètre 


Quand on considère le déplacement du centre d’inté- 
rêt de ces mesures d'indices de réfraction, on voit donc 
que, d’abord attaché à des zones de grande transpa- 
rence des substances, pour des fins essentiellement pra- 
tiques (construction de spectrographes), pour des 
raisons théoriques et de nouvelles applications, il 
englobe maintenant les régions d'absorption. De nou- 
velles techniques ont été donc en même temps mises 
au point pour s'adapter au cas des corps absorbants. 
Ce sont ces méthodes, tant anciennes que nouvelles, 
que nous allons maintenant passer en revue, en nous 
hmitant, dans chaque cas, aux résultats récents. 


I. Méthode du prisme. — Chaque fois que l’on a 
affaire à un corps transparent, et se présentant sous 
forme de cristaux de grandes dimensions, c’est la 
méthode du prisme qui se montre la plus précise, soit 
que l’on opère au minimum de déviation, soit que 
l’on utilise l’autocollimation. Ces procédés anciens, 
simples et bien connus, ont été l’objet d’amélio- 
rations de détail mais ils conduisent à une précision 
accrue, S'il s’agit seulement d'étendre au proche 
infrarouge les méthodes classiques dans le visible, il 
suffit de mettre à la sortie de l’oculaire un conver- 
tisseur d'image [11] ; on arrive ainsi à 1,1 & ou même 
à 1,4 &. Mais nous nous intéresserons ici surtout aux 
méthodes spéciales à l’infrarouge et permettant d’aller 
beaucoup plus loin. 

Rodney et ses collaborateurs transposent sim- 
plement dans l’infrarouge la méthode classique dans le 
visible, du prisme monté sur la plateforme d’un gonio- 
placé au minimum de déviation ; à l’extré- 
mité de l’oculaire, l’œil se remplace par le récepteur 
infrarouge ; par suite du temps assez long que deman- 
derait la recherche directe du minimum de déviation 
pour chaque radiation, on simplifie les mesures en 
remarquant que, en passant d’une radiation à la sui- 
vante, la variation de la déviation représente le double 
de celle de l’angle d'incidence sur le prisme : les mou- 
vements des deux bras du goniomètre se lient donc 
l’un à l’autre par des engrenages [61] ou par d’autres 
systèmes mécaniques plus anciens. 

On a mesuré ainsi l'indice du bromure de césium [611, 
de la silice fondue [62], du bromoiodure de thal- 
hum [63], du saphir synthétique [45] et du trisulfure 
d’arsenic [64]. La précision atteint quelques unités 
de la 5° décimale. 

La même méthode, mais avec autocollimation, a été 
utilisée à la mesure de l’indice de quelques verres dans 
le proche infrarouge [17], [7] et de sulfures [14] dans 
leur zone de transmission. C’est par autocollimation 
également, en augmentant considérablement la préci- 
sion qu'ont été mesurés les indices de 193 verres 
d'optique entre le visible et 1,7 p [55]. L'erreur ne 
dépasse pas 6.10—5 sur l'indice. 

Toutes les méthodes précédentes supposent un mon- 
tage particulier et permanent. Quand on a à effectuer 
un petit nombre seulement de mesures, il peut être 
avantageux d'utiliser un spectrographe ordinaire en 
lui apportant quelques modifications. Dans cet ordre 
d'idées, nous relevons trois méthodes : MacAlister, 
Villa et Salzherg [3] installent le prisme à étudier à la 
place du prisme ordinaire dans un spectrographe infra- 
rouge Perkin-Elmer et l’étudient par autocollimation. 
Les déterminations se font rapidement avec une 
précision encore bonne (2 unités de la 4e décimale) ; 
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cette méthode a été appliquée à l’étude de verres 
d'optique [12], de la silice fondue à haute tempé- 
rature [79], du sélénium et du germanium [65], [66] ; 
ces deux derniers corps présentent des difficultés sup- 
plémentaires car, étant opaques dans le visible, aucun 
réglage visuel ne s’avère possible, et la méthode 
d’autocollimation reste précieuse dans ce cas ; leur 
dispersion à peu près nulle dans certains domaines 
oblige d’autre part à utiliser un dispositif dispersif 
supplémentaire. D’autres auteurs recourent au même 
spectrographe, mais en remplaçant la boîte contenant 
la source par un montage goniométrique comportant 
le prisme étudié au minimum de déviation [18]. 
L'indice du sulfure de cadmium entre 1 et 15 1, a été 
déterminé de cette manière. 

Nous mêmes avons réalisé, au Laboratoire de 
Recherches Physiques, un réfractomètre pouvant 
s’adjoindre à un spectrographe Perkin-Elmer simple 
faisceau, en conduisant à la sortie de celui-ci, les ra- 
diations monochromatiques dans un montage de 
Wadsworth avec prisme à étudier au minimum de 
déviation et en faisant rentrer ensuite ces radiations à 
l’intérieur du spectrographe pour les détecter. Nous 
avons réalisé, avec cet appareil, des mesures d’indice 
de verres de la maison Parra Mantois [87]. 

La mesure des indices du tellure [25], [8] jusqu’à 
14 x et de l’antimoniure de gallium [16] se réfère aussi 
à la méthode du prisme, avec des variantes. 

Enfin la dispersion de quelques gaz comme HCI1[32], 
HBr et CF, [67], [68], se mesure par la méthode du 
prisme creux, déjà utilisée autrefois par Rollefson. 
Dans ce dernier cas, l’intérêt des mesures ne se relie 
plus à lPutilisation pratique des substances étudiées, 
comme dans les cas qui précèdent, mais il s’agit alors 
de calculer l'intensité de bandes d’absorption, extrè- 
mement fortes, et d’en déduire ensuite, avec l’aide 
d’une hypothèse sur la forme de la fonction potentiel, 
les valeurs du moment électrique de la molécule et de 
sa dérivée première par rapport à la distance inter- 
nucléaire. 


IT. Méthode de la réflexion totale. — Introduite 
dans l’infrarouge par Seegert, elle avait été utilisée 
par l’un d’entre nous à la mesure de l’indice de quelques 
liquides entre le visible et 3 u. Nous venons de la 
reprendre depuis peu, d’abord avec un réfractomètre 
constitué par des prismes en verre, ensuite par des 
prisme en pentaseléniure d’arsenic, avec lesquels nous 
pourrons effectuer des mesures jusqu’à 18 4 [87]. 

Dans un travail récent, Jaffé et ses collabo- 
rateurs [31], choisissent à nouveau la méthode en lui 
apportant quelques améliorations, en particulier en. 
étendant le domaine étudié jusqu’à 6 , grâce à lutili- 
sation de deux prismes en trisulfure d’arsenic ; ils 
déterminent ainsi l’indice de quelques liquides orga- 
niques. La méthode se montre précieuse, en permettant 
des mesures d'indices relativement précises 

(An =510) 
dans des régions de faible absorption, où la méthode 
du prisme creux ne peut plus s'appliquer. 

Il semblait qu'elle cessait d’être utilisable lorsque 
l'indice d’absorption x dépassait la valeur 0,01, le 
phénomène de réflexion totale n’existant plus à pro- 
prement parler. Malgré ce grave inconvénient, Jaffé 
dans un article très récent [33], utilise tout de même 
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la mesure de la lumière transmise à l'étude de la varia- 
tion de l'indice, à travers des bandes d'absorption pas 
trop fortes de certains liquides. Pour cela, il réduit 
l'épaisseur de la couche liquide à des valeurs très 
faibles, de l’ordre de grandeur de la longueur d’onde, 
et il introduit des caleuls correctifs, pour tenir compte 
de la valeur de l'indice x. La précision devient naturel- 
lement moindre (An — 0,03 pour x — 0,03) et néces- 
site, pour les corrections, Putilisation d’un calculateur 
électronique. 


III. Méthodes interférentielles. — Dans l’infra- 
rouge, les substances présentent, pour la plupart, de 
nombreuses bandes d'absorption, et les méthodes pré- 
cédentes, les plus précises, correspondent à un champ 
d'application restreint. Dès que l’on se trouve au voisi- 
nage des bandes d'absorption, elles se montrent inuti- 
lisables et, dans ces conditions, les méthodes interfé- 
rentielles, moins précises en général, deviennent inté- 
ressantes ; il en est de même dans le cas de substances 
solides transparentes, mais de très haut indice. 

L'étude des spectres cannelés, appliquée depuis 
longtemps, a été utilisée dans plusieurs travaux récents 
avec de nombreuses améliorations, tant au point de 
vue pratique que dans les calculs. 

S'il s’agit de solides, la méthode consiste à tailler 
une lame mince à faces parallèles, à mesurer son épais- 
seur, et à pointer les maximums des franges d’inter- 
férences obtenues. La formule fondamentale 
2ne — KA (1) valable tant que la substance reste 
transparente ou légèrement absorbante, permet alors 
de calculer n, connaissant K. En réalité, il faut bien 
souligner que la détermination exacte de K ne se 
montre en général pas possible avec un seul spectre, 
sauf si l’indice nr est déjà connu en un point de l’inter- 
valle étudié : Moss et ses collaborateurs [48] insistent 
sur cette question que certains auteurs semblent négli- 
ver [52] et nous-mêmes avons montré, dans le cas des 
liquides [83], qu’il était nécessaire de comparer plu- 
sieurs spectres, correspondant à des épaisseurs diffé- 
rentes. 


Par cette méthode ont été déterminés les indices de : 


l’antimoniure d’indium [48] entre 7 et 20 , du sulfure 
de zinc entre le visible et 15 1 [51], du germanium [52], 
dans le proche infrarouge, ces dernières mesures ayant 
suscité les réserves citées plus haut, de l’oxyde de 
titane [13] et du sélénium entre 9 et 23 u [8]. 

S'il est diflicile de tailler des lames de l’épaisseur 
désirée, il devient quelquefois possible, par contre, de 
déposer le corps en couche mince par évaporation ther- 
mique, en prenant toutes précautions pour aboutir à 
une couche propre et non lacunaire : l'indice du sélé- 
nium à été mesuré ainsi dans le proche infrarouge [41]. 

Nous nous étendrons un peu plus longuement, à titre 
d'exemple, sur le cas des liquides, que nous venons 
d'étudier personnellement [81, 83, 84] et qui se montre 
plus délicat que celui des solides. 

Le premier problème, qui se pose, concerne l’interfé- 
romètre et c’est là qu’ont été réalisées les améliorations 
les plus importantes, correspondant à l’utilisation dans 
l’infrarouge de nouveaux corps transperents de haut 
indice et au dépôt, par évaporation thermique, d’un 
nombre de plus en plus grand de diélectriques variés. 
Kagarise [35] utilise un interféromètre, constitué par 
des lames de germanium et nous-mêmes [84] une cuve 
à épaisseur variable avec des faces de bromoiodure de 


No 


thallium. On peut aussi prendre des lames d’une sub- 
stance de bas indice et les recouvrir d’une couche 
mince d’un diélectrique transparent de haut indice 
comme le pentaséléniure d’arsenic par exemple [81]. 
Dans ce dernier cas, l'épaisseur de ces couches n’est pas 
indifférente, et il faut prendre garde à certains phéno- 
mènes d’interférences additionnels [82]. 

Au point de vue des calculs, et comme dans le cas 
des solides, le problème le plus diflicile à résoudre 
reste la détermination de l’ordre d’interférences K 
dans la formule (1). L'étude d’un seul spectre cannelé 
ne nous fournit qu'une relation différentielle entre la 
variation de X et celle de n. La détermination absolue 
de À exige le dépouillemènt de plusieurs spectres 
cannelés, obtenus en faisant varier une des grandeurs 
de la relation (1). Dans cet ordre d’idées deux nouveaux 
procédés ont été proposés : Kagarise [37] opère à épais- 
seur constante et modifie » d’une façon continue, en 
partant d’un liquide d'indice connu, le sulfure de 
carbone, par exemple, et mélangeant à celui-ci le 
liquide d’indice cherché, èn proportion de plus en plus 
grande, jusqu’à obtenir ce dernier à l’état pur. On 
suit alors une cannelure d’ordre donné par continuité 
d’un spectre au suivant. Nous-mêmes [83], opérant au 
contraire toujours avec le liquide pur, faisons varier 
l'épaisseur d’une façon continue, grâce à l’utilisation, 
comme interféromètre, d’une cuve à épaisseur variable. 
L'une ou l’autre de ces méthodes permet de déter- 
miner À sans ambiguité. 

Malgré ces perfectionnements, la méthode interfé- 
rentielle des spectres cannelés fournit une précision 
bien moindre que les précédentes. L’une des raisons 
provient de l’erreur sur la mesure de l’épaisseur e, 
l’autre de ce que l’on travaille sur un spectre continu ; 
le repérage des longueurs d’onde est alors beaucoup 
moins précis que si l’on opère avec des raies étalons. 
(Rank et ses collaborateurs [52] ont toutefois obvié à 
cette dernière difficulté dans le proche infrarouge en 
enregistrant en même temps spectre cannelé et raies 
d'émission.) La précision sur n ne dépasse pas, en 
général, quelques unités de la 3e décimale. 

Les mesures d'indices de Kagarise et ses collabora- 
teurs [36] ont été effectuées dans le but de déterminer, 
par ce moyen détourné mais précis, les intensités de 
bandes très fortes de liquides. 

Enfin le procédé interférentiel a été appliqué avec 
succès à la mesure d'indices de gaz. Rank et ses colla- 
borateurs [53, 54] ont déterminé ainsi, avec une haute 
précision, l'indice de l’air dans le proche infrarouge, 
Legay et collaborateurs [43, 48], l'indice du gaz 
chlorhydrique et du gaz carbonique à l’intérieur même 
des bandes d'absorption, entre les raies de vibration- 
rotation. 

La méthode interférentielle s’étend sans difficulté au 
cas où la substance est légèrement absorbante. Tant 
que x reste petit par rapport à n, la position des 
maximums et des minimums demeure inchangée et 
correspond à des ordres d’interférences entiers. Par 
contre la transmission se trouve diminuée. La position 
des franges permettra donc encore de calculer n, comme 
précédemment, et la mesure de la transmission maxi- 
mum, par exemple, conduira au calcul de x. Genzel et 
ses collaborateurs [20, 21] ont déterminé ainsi les cons- 
tantes optiques de LiF, NaCI, KCI et KBr dans l’infra- 
rouge lointain, au delà de leur bande d’absorption. 
Hs trouvent, comme il fallait s’y attendre, une valeur 
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de 7 constante et égale à la racine carrée de la cons- 
tante diélectrique statique. Markov et Lindstrem [46] 
appliquent la même méthode au bismuth, Huldt et 
Staflin [30] au sulfure de zinc et au germanium en 
couches très minces ; ces derniers montrent que les 
constantes optiques varient, dans une certaine mesure 
avec les conditions d’évaporation et que ces évapora- 
tions doivent donc être effectuées avec beaucoup de 
SOIN. 

Enfin, à ces méthodes interférentielles, on peut 
rattacher des mesures, non pas d’indices mais de 
biréfringences, au moyen de spectres cannelés de biré- 
fringence [80, 69], en lumière polarisée, 


IV. Méthodes par réflexion. — Nous avons vu que 
l'intérêt porté aux mesures d'indices de réfraction, 
d’abord réservé aux substances transparentes, s'était 
étendu ensuite aux zones d'absorption de ces mêmes 
substances et même aux corps complètement absor- 
bants comme les métaux. Sauf cas très particuliers, 
aucune des méthodes précédentes ne convient dans ces 
conditions, puisqu'elles impliquent toutes la traversée 
de la substance par le rayonnement, et qu’on ne dis- 
pose alors, à la sortie, que d’une énergie trop faible 
pour effectuer des mesures. L’énergie réfléchie par 
contre devient en général considérable, et presque 
toutes les mesures d'indices de corps très absorbants 
sont fondées sur des études par réflexion. 

Une difficulté se pose dès l’abord : toute propriété 
optique de la substance est fonction, non seulement de 
l’indice de réfraction n mais aussi de l’indice d’absorp- 
tion x. Lorsque x reste petit par rapport à n, inférieur 
à 0,1 par exemple, ce qui peut déjà correspondre à une 
absorption assez forte, le pouvoir réflecteur ne dépend, 
en réalité, que de x à l’approximation obtenue en 
général, et la mesure permettra de calculer n direc- 
tement. Par contre, lorsque xn’apparaît plus petit par 
rapport à n,les contributions de n et x deviennent 
impossibles à séparer dans la réflexion, et nous aurons 
deux grandeurs à déterminer au lieu d’une seule. Deux 
manières différentes de considérer le problème s'offrent 
alors. Nous pouvons considérer n et x comme deux 
grandeurs quelconques et effectuer deux mesures indé- 
pendantes qui nous fournissent deux relations numé- 
riques entre 7» et x ; mais nous pouvons aussi, partant 
d’une étude théorique des phénomènes d'absorption 
et de dispersion, chercher, à priori, une relation entre n 
et x et n’effectuer ensuite qu’une seule mesure pour 
déterminer ces deux quantités. (Il est important de 
remarquer ici que les deux mesures ou la mesure 
unique s'entendent à longueur d’onde déterminée.) 

Nous allons maintenant passer en revue les différents 
procédés utilisés pratiquement en nous y étendant un 
peu plus largement que sur les précédentes ; ces 
mesures par réflexion sont, en effet, beaucoup moins 
classiques que celles par réfraction, ou même par inter- 
férences et elles ont pris ces dernières années une grande 


extension. 


A. SUBSTANGES MOYENNEMENT ABSORBANTES. — Le 
pouvoir réflecteur s'exprime par la formule 
(n — 1)? + x? 
GAP X? 


Étant donné que la précision sur /? reste en cénéral, 
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si x est inférieur à 0,1. Dans ce cas, x pourra être 
mesuré par ailleurs, par transmission, par exemple. 

La méthode se recommande par sa simplicité, mais 
non par sa précision très médiocre ; elle ne s’utilisera 
que pour arriver rapidement à une valeur approchée 
d’un indice [15], ou lorsque les méthodes précédentes 
deviennent inutilisables. Caldwell et Yan [8] ont déter- 
miné, de cette manière, lindice du sélénium entre 30 
et 150 x. Les anciennes mesures de Centano sur l’eau 
reposent sur ce principe. Ces derniers temps, la méthode 
a été reprise par M. Caméo [10] pour des liquides, en 
dehors de leurs bandes d'absorption forte. Enfin des 
composés organiques, se présentant sous forme de 
poudre, peuvent s’étudier d’une manière semblable, 
en les faisant fondre puis solidifier au contact d’une 
lame d'indice connu [9]. + 

Les constantes optiques du chloroforme se trouvent 
dans le travail de Lisitsa et Tsiashchenko [44] ; x se 
mesure par absorption et n par réflexion mais aucun 
détail n’est donné sur la méthode utilisée par réflexion, 
pour éliminer Pabsorption de la vapeur. 4 


voisine de 1 %,, la mesure de ? donnera directement n, 


B. SUBSTANCES TRÈS ABSORBANTES. — 10 Méthodes 
utilisant deux mesures à chaque longueur d'onde. — Sans 
faire aucune hypothèse sur les relations qui peuvent 
exister entre n et x, on effectue deux mesures indé- 
pendantes de grandeurs liées à n et x. Plusieurs mé- 
thodes ont été proposées récemment : 

a) Une des techniques les plus simples consiste à 
mesurer la transmission en couche très mince et le 
pouvoir réflecteur sous imcidence normale. Ce procédé, 
qui se rapproche de celui indiqué plus haut en A, a été 
rendu possible dans le cas des substances très absor- 
bantes par le fait que l’on sait réaliser maintenant, 
pour de nombreux corps, par évaporation thermique 
sous vide, des couches aussi minces, donc aussi peu 
absorbantes que l’on désire. Mais les propriétés 
optiques des couches très minces se montrent souvent 
différentes de celles du même corps à l’état massif, et 
il faut être très prudent dans la généralisation des 
résultats obtenus par de telles méthodes. Dans cet 
ordre d’idées, on peut citer les travaux de Hall et 
Ferguson [24] sur les sulfures de zinc et de cadmium, 
de Yoshinaga et Oetjen [88] sur l’antimoniure d’indium 
et de Hass [27] et collaborateurs sur les oxydes et 
fluorures de terres rares. 

b) Si l’on utilise uniquement des mesures par 
réflexion, il faut exécuter deux mesures dans des condi- 
tions différentes, Simon [73, 74], par exemple, utilise 
deux angles d'incidence, 209 et 700 ; le montage appa- 
raît assez délicat à réaliser et les calculs numériques 
sont fastidieux ; on peut toutefois les simplifier en 
utilisant toujours les mêmes angles d'incidence et 
construisant alors des abaques une fois pour toutes. 
D'une manière générale, les déterminations par 
réflexion ayant une précision médiocre, dépassant rare- 
ment 1 %, permettent souvent l’utilisation de gra- 
phiques ou abaques qui simplifient beaucoup les 
calculs. Cette méthode présente toutefois un incon- 
vénient plus grave, dans le cas de substances aniso- 
tropes, l'interprétation des résultats, pour l'incidence 
de 709, devenant alors très pénible. La méthode de 
Simon a été utilisée ensuite par d’autres auteurs : 
Simon, Mac Mahon [75] et Reitzel [56] l’emploient 
pour localiser avec précision la position des bandes 
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d'absorption du quartz et de différentes formes de la 
silice cristallisée ou fondue et étudier le passage d’une 
forme à une autre, quand la température varie. 
D’autres auteurs [5] mesurent les constantes optiques 
de l’arséniure de gallium, et en déduisent la masse 
effective des électrons de conduction. Enfin 
Neuroth [49] et Heilman [29] étudient, de la même 
manière, des verres entre 3 et 12 y et le fluorure de 
lithium entre 15 et 36 y. 

Avery [4], d'autre part, préfère une méthode ana- 
logue, mais simplifiée, en se contentant de mesurer des 
rapports d'énergie réfléchie dans deux plans de polarisa- 
tion différents. Roberts [57] détermine les constantes 
optiques du nickel et du tungstène en mesurant, sous 
grande incidence, les énergies réfléchies dans deux états 
de polarisation différents et en tire une étude théorique 
intéressante des électrons de conduction dans ces 
métaux. L’étude des métaux est rendue difficile par 
le fait que leur pouvoir réflecteur très élevé, voisin 
de 1, exige, dans sa détermination, une précision supé- 


rieure: à ce que l’on obtient habituellement. Des” 


auteurs [ 90, 71, 72] obvient à cette difficulté grâce à 
des réflexions multiples et déterminent, de cette ma- 
nière, les constantes optiques de nombreux métaux à 
très haut pouvoir réflecteur. 

c) Pour éviter les difficultés expérimentales et les 
calculs pénibles qui résultent de l'utilisation de grands 
angles d'incidence, surtout quand il s’agit de corps 
anisotropes, plusieurs méthodes à incidence normale 
ont été proposées à peu près simultanément. Abelès [1] 
avait indiqué la possibilité, aans l’étude des métaux, 
de recouvrir ceux-ci d’une couche mince transparente, 
à la fois pour les protéger de loxydation superficielle 
et pour faire varier le pouvoir réflecteur. Nous avons 
proposé et utilisé de notre côté [85, 86] une méthode 
fondée sur la détermination du pouvoir réflecteur, sous 
incidence normale, d’une part de l'échantillon solide 
étudié, et, d’autre nert, du rôme échantillon recou- 
vert, par évaporation, d’une couche mince d’un corps 
d'indice connu. Le dispositif expérimental apparaît 

.ainsi très simple et l'interprétation des résultats 
conduit à des calculs élémentaires. S'il s’agit d’un 
liquide, une variante de la méthode [10] consiste à 
interposer une lame épaisse d’un corps d'indice connu. 

La méthode de Kravet, fondée sur un principe un 
peu analogue, à été utilisée pour des métaux plongés 
dans des liquides d’indices variés et connus [76]. 

Enfin, Kendall et Turner [38] semblent, malgré le 
peu de détails de leur note, avoir effectué sur le fer, 
en couche mince, un travail tout à fait semblable aux 
précédents. 

20 Méthodes utilisant une seule mesure du pouvoir 
réflecteur à chaque longueur d'onde. — a) L'indice de 
réfraction et l’indice d'absorption ne constituent pas 
deux grandeurs indépendantes ; c’est ce que l’on sous- 
entend implicitement en les faisant entrer dans un. 
indice complexe N — n — jx. On suppose ainsi que N 
est une fonction analytique de la fréquence. Ceci 
revient, en fait, à supposer qu'il existe une relation 
linéaire entre le champ électrique et le déplacement 
électrique à l’intérieur de la substance considérée, Les 
conclusions, que l’on peut tirer de cette hypothèse, 
ont été étudiées ces dernières années, d’abord au point 
de vue électrique [6, 19], ensuite au point de vue 
optique [47]. 

S'il s’agit de l'étude du pouvoir réflecteur, qui nous 
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intéresse ici, nous savons que le rapport R& des ampli- 
tudes des champs électriques réfléchi et incident est, 


à chaque instant, un nombre complexe, si la substance 

N —1 : 
onsidérée est absorbante : R = 57 — r ef®. 
considérée est abs ù Ni 


Le pouvoir réflecteur mesuré s’écrit : À = r?. 

Supposant que R est une fonction analytique de la 
fréquence, on trouve qu'il existe entre r et @ une rela- 
tion intégrale. Autrement dit, si l’on connaît la varia- 
tion de À pour toutes les fréquences, on peut en dé- 
duire, par le caleul numérique d’une intégrale, le dépha- 
sege o en chaque point. 

Pratiquement certains domaines de fréquences pos- 
sèdent une influence prépondérante par rapport à 
d’autres, et on peut se limiter à un intervalle fini. 

Cette méthode, très intéressante théoriquement, 
permet de déterminer » et x à partir de la seule mesure 
du pouvoir réflecteur À dans un certain domaine de 
longueurs d'ondes. Elle a été appliquée, dans le cas 
de cristaux, par Robinson et Price [58, 59, 60], Hass 
et Sutherland [26], Abelès et Mathieu [2]. Elle se 
montre particulièrement précieuse dans le cas des 
cristaux ioniques, qui présentent un très petit nombre 
de bandes. Son inconvénient réside dans la longueur 
du calcul numérique des intégrales définies qu’eile 
implique ; l'emploi d’un calculateur électronique sim- 
plifie alors considérablement ces calculs. 

Dans tous les travaux précédents, il s'agissait de 
déterminer la position et la forme exactes des bandes 
d'absorption. Dans le cas des cristaux, en effet, la 
détermination directe en reste souvent impossible, à 
cause de la grande intensité des bandes et de leur 
variation avec la direction du cristal. La détermination 
simultanée de n permet ensuite de comparer les valeurs 
obtenues à celles que donne la théorie de Drude, par 
exemple. La même méthode a été appliquée par 
Gottlieb [22], par Philipp et Taft [50] à la bande 
d'absorption du germanium, et par Kislovskni [39, 40] 
à des liquides et des cristaux. 

b) Inversement, partant d’une théorie comme celle 
de Drude, qui nous fournit une relation entre n et x 
on peut déterminer facilement nr et x à partir de la 
seule mesure du pouvoir réflecteur, et vérifier ensuite 
si les hypothèses faites sont exactes, aux erreurs expé- 
rimentales près. Cette méthode, moins générale que la 
précédente, puisqu'elle implique une théorie atomique 
des propriétés optiques, se montre d’une exécution 
beaucoup plus simple. De nombreux auteurs l’ont 
appliquée avec succès, C. Haas [23] à la calcite et au 
nitrate de sodium, Spitzer et collaborateurs [77, 78] 
aux différentes variétés de carbure de silicium. Si l’on 
veut ajuster les constantes rentrant dans les formules 
de dispersion, de manière à obtenir la meilleure concor- 
dance possible avec les résultats expérimentaux, il faut 
toutefois se lancer dans de longs calculs numériques de 
tâtonnement. 

Encore beaucoup plus simple est la méthode essayée 
par Simon [73], et qui consiste, parlant le langage 
électrique, à appliquer les résultats de la transmission 
du champ électromagnétique le long des fils, et à 
utiliser les graphiques de Smith. Cette méthode, extrê- 
mement simple, part de l'hypothèse de vibrateurs 
harmoniques, et suppose que l'indice complexe du 
milieu absorbant varie avec la fréquence de la même 
manière que l’impédance complexe d’un circuit réso- 
nant ; elle conduit toutefois à des résultats SrOSSiè- 


2. 
Nout 


rement approchés, et ne peut être considérée comme 
une véritable méthode de mesure des constantes 
optiques. 


V. Filtres de Christiansen. — Nous ne reviendrons 
pas sur le principe de la méthode des filtres de Chris- 
tiansen déjà exposé précédemment [421]. La question 
vient d’être reprise récemment [34, 87] et nous en 
tirerons seulement quelques conclusions, relatives à la 
mesure des indices par ce procédé. 

Comme on le sait, si l’on met en suspension une 
poudre dans un liquide, ou dans un gaz, on observe 
un maximum de transmission à la longueur d’onde 
pour laquelle les indices des deux corps se présentent 
comme égaux ; si l'indice du liquide est connu ou s’il 
s’agit d’air dans les conditions normales, on en déduit 
immédiatement l'indice du solide, et inversement. Tou- 
tefois la précision de cette détermination se montre 
très faible et atteint 1 % dans les meilleures condi- 
tions. En effet, si le liquide et la poudre sont parfaite- 
ment transparents dans la région considérée, les deux 
courbes de dispersion possèdent des pentes voisines 
et, dans ces conditions, le filtre obtenu apparaît très 
large avec un maximum difficile à pointer. Au con- 
 traire, si l’on se trouve au voisinage d’une bande 
d'absorption de l’un des corps, les courbes de disper- 
sion peuvent présenter des pentes très différentes et 
le filtre se montrer très étroit ; par contre l’absorption 
propre de la substance risque alors de déplacer le maxi- 
mum de transmission vers l’extérieur de la bande, et 
de donner une erreur systématique dans l'indice. Dans 
ce dernier cas, l’utilisation de plusieurs filtres d’épais- 
seurs différentes permet d'observer le phénomène et 
d’y remédier dans une certaine mesure. 

La méthode des filtres de Christiansen a tout de 
même été utilisée dans certains cas particuliers, en 
complément d’autres méthodes plus précises : dans 
l’étude des spectres cannelés par exemple, dans les 
domaines où l’on ne possède aucune idée, à priori, de 
la valeur de l’ordre d’interférence Æ, la réalisation 
d’un filtre de Christiansen fournit alors une valeur 
approchée de l'indice, qui peut être précieuse dans le 
choix de l’ordre d’interférences. S'il s’agit de filtres de 
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poudres dans Pair, on saura ainsi si l’indice du solide 
passe ou non par la valeur 1 [87]. 


Conelusion. — Noûs avons vu, dans ce qui précède, 
le développement considérable, pris par les mesures 
d'indices dans l’infrarouge ces dernières années. De 
nombreuses méthodes ont été proposées, applicables à 
des corps aussi bien absorbants que transparents, et 
dans des états physiques variés. 

Tant qu’il s’agit de substances solides, transparentes 
dans un domaine étendu de longueurs d’ondes, la 
vieille méthode du prisme s'avère toujours la meil- 
leure; s’il y a absorption faible, les méthodes interfé- 
rentielles deviennent utiles et enfin, si l’on s’adresse à 
des corps très absorbants, c’est aux procédés par ré- 
flexion qu’il faut recourir. Le cas des liquides se pré- 
sente comme plus délicat, car ils possèdent tous des 
bandes d’absorption plus ou moins intenses dans 
l’infrarouge, et les méthodes de réflexion totale, de 
spectres cannelés et de réflexion sous incidence normale 
présentent toutes des avantages et des inconvénients ; 
on rencontre encore quelques cas où aucun de ces 
procédés ne semble vraiment satisfaisant. Enfin, s’il 
s’agit de déterminer l’indice de gaz, la méthode interfé- 
rentielle et celle du prisme creux donnent toutes deux 
d'excellents résultats. 

On peut donc dire que chaque cas particulier a 
reçu, à quelques exceptions près, une solution appro- 
priée. Il reste, pour les appliquer, un champ encore 
vaste, aussi bien en ce qui concerne le domaine spectral 
que la variété des substances à étudier. Ce sera sans 
doute la tâche des prochaines années. On constate, en 
effet, actuellement une disproportion vraiment trop 
grande entre le nombre de spectres d'absorption (ou de 
réflexion) infrarouge et la connaissance des indices de 
réfraction dans le même domaine. On explique aisé- 
ment d’ailleurs ce résultat par suite de la facilité des 
mesures dans le premier cas et de leur difficulté dans 
le deuxième, mais nous avons essayé de souligner tout 
l’intérêt que présente la détermination de la dispersion 
des substances et les possibilités variées que l’on pos- 
sède pour aborder cet important problème. 


Manuscrit reçu le 22 juin 1960. 
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LA RÉDACTION 


POLARISATION NUCLÉAIRE 
PAR « DOUBLE EFFET », 
EN PRÉSENCE DE DIFFUSION, 
DANS LE CAS 
OU LA RAIE DE RÉSONANCE ÉLECTRONIQUE 
EST INHOMOGÈNE 


Par Jean-Loup MOoTcHANE et Jean UEBERSFELD, 


École Supérieure de Physique et Chimie de Paris, 
Faculté des Sciences de Besançon. 


Dans une communication antérieure [1] on a montré 
que la polarisation par « double effet » des noyaux 
d’un solide contenant des impuretés paramagnétiques 
était affectée par le caractère « inhomogène » de la 
raie de résonance des spins électroniques. 

Le coefficient d'augmentation de la polarisation 

nucléaire p correspondant à l’une des 2 transitions de 
fréquence v, + vx, où v, représente la fréquence cen- 
trale de la fonction de répartition G(v) des fréquences 
de résonance électronique ve, et vx la fréquence de réso- 
nance nucléaire dans un champ magnétique cons- 
tant Æ,, est donné par l’expression : 


Tz + À Te Vo 


=, L'as)-12 D—1{s) (1 
ERA 1 1, ni - as) (s) (1) 
avec 
# nn D 
Ds) ol 1 PR = { RFA as) 
3B HE NE ie B (2) aies \ 
+) 4 +00 + (0) (+ 05) 
et 
A œ 
1 4 L % 
ho D A | 
On appelle : 


L,, la polarisation dynamique moyenne des noyaux, 
15, la polarisation nucléaire à l'équilibre thermo- 
dynamique, 
27 T,, l'inverse de la demi-largeur à mi-hauteur de 
la fonction de répartition G(v). | 
t, et t, sont les temps de relaxation électronique, 
spin-réseau et spin-spin, ye le rapport gyroma- 
gnétique de l’électron, À, l’intensité du champ 
magnétique tournant de fréquence v = v5 = vx. 


w = 7, et w. sont les probabilités de relaxation pure- 


ment électronique et nucléaire. 
Woo est la probabilité de relaxation aux fréquences 
Ve Æ YN. 
Il est possible d'étendre les calculs antérieurs [1] au 
cas où l’on étudie la polarisation des noyaux d’un 
fluide au contact d’un solide paramagnétique poreux. 


Le modèle moléculaire utilisé est le suivant [2]. On 
considère un pore du solide, constitué par exemple, par 
une sphère creuse de volume V et de rayon À qui 
contient le fluide. On suppose que les centres para- 
magnétiques sont situés à la surface de la sphère et 
qu’il existe au voisinage de cette surface une couche de 
fluide adsorbé. 

On peut classer les centres paramagnétiques en 
2 groupes. Le premier contiendra les centres para- 
magnétiques dont les fréquences de résonance w sont 
telles que ve + ve, ve — vx ou w soit égal à la fré- 
quence y de la transition induite par le champ magné- 
tique tournant d'intensité A.. 

Ce sont ces centres qui sont responsables de la pola- 
risation nucléaire. On appellera Vr le volume occupé 
par les noyaux de la couche du fluide adsorbé qui sont 
liés par une interaction dipolaire statique à ces centres 
paramagnétiques. 

Le deuxième groupe renferme les autres centres para- 
magnétiques. Ils ne participent que faiblement à la 
polarisation nucléaire mais constituent une cause de 
relaxation importante pour les noyaux en interaction 


: avec eux. 


On appellera Vs le volume occupé par les noyaux 
de la couche du fluide adsorbé qui ne sont pas compris 
dans Vo. 

On considérera 3 régions dans le fluide : 

— Les régions 1 et 2 qui compréndront respecti- 
vement les noyaux occupant les volumes Vr et Vs. 

— La région 3 dans laquelle les noyaux, par suite 
de leur mouvement rapide, ne sont soumis à aucune 
interaction statique par les centres paramagnétiques. 

Les molécules du fluide passent constamment de la 
région 3, aux régions {À et 2 et réciproquement. 

On désignera respectivement par 7, et + (temps 
d'accrochage) les temps pendant lesquels une molécule 
venant de la région 3 reste dans l’une des régions { et 2. 
On voit que pour que l’on puisse parler d’une inter- 
action dipolaire statique entre les spins électroniques 
et nucléaires, dans les régions 1 et 2 il faut que les 
temps 7, et +, soient plus longs que la durée de vie des 
états quantiques intervenant dans les transitions de 
fréquence v. 

Le temps de relaxation spin réseau 7, des noyaux 
dans l’ensemble du fluide est donné par la formule (2) : 


RS EE" 
PAT ep 


T;» désigne le temps de relaxation spin-réseau des 
noyaux dans la région f. 

Le temps de relaxation spin-réseau T°, des noyaux 
dans la région 3 compte tenu de la relaxation dans la 
région 2, est donné par la formule de Torrey [2] 


il il Vs il 


FR er CN PT PER 


où Ty et T's sont respectivement les temps de rela- 
xation spin-réseau des noyaux dans les régions 2 et 5, 
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L'équation d'évolution en fonction du temps de la 
polarisation nucléaire /;, en présence de diffusion molé- 
culaire (Équation de Torrey et de Bloch [3]) s'écrit 
dans la région à. 

PR ae € 
dt Ti 


où D est la constante de diffusion. 

Si la diffusion est assez rapide, c’est-à-dire si 
(DT'}2 5 R, la polarisation nucléaire Z, est une cons- 
tante en régime permanent. 

IL est relié à Z, polarisation nucléaire moyenne dans 
la région 1, par la condition de conservation de la pola- 
risation en régime permanent, à la frontière des 
régions À et 3 : 


I, dans le cas où la raie de résonance électronique est 
inhomogène est donné par lexpression : 


7 C2) ’ / 2 Q 
Ba = [7 Tale) Ge) de (3) 


où Z,(v{) est la polarisation des noyaux de la région 1 
en interaction avec les spins électroniques dont la fré- 
quence de résonance est ve. x 

G(ve) est la fonction de répartition des fréquences vw. 


I, calculé par des méthodes exposées antérieu- 
rement [4] est donné par la relation : 
HS SOUDE ARS Sir 
a ee E E 2(oo + Y) Jalve) + ALOT (4) 
avec 
/ 2 l ñ 
vel D’ (oi — Bi) (2 + 200) 
7 4 1 ER) l R 95 
fa (ve) = D’ (2 + 249 + &i + Pi 20;) 
D' = (20 L 2y + =) 12 + 240 + œ1 + Pi + 20:) 
0 WT D 0 1 1 1 
+ (oi + Ba) (2 + 200 + 281) + kon Br 


S, représente la polarisation électronique à l’équilibre 
thermodynamique. 

Wos We, eb ws, probabilités de transition induites 
par le champ magnétique tournant /7;, aux fréquences 
respectives v6 + vw, ve — vx et w sont données par les 
relations : 


Gi — Ep SÉ(Ve, Y — YN) 
A 
ECS AN € 


Î(ve, v) est le profil de la courbe de résonance d’un 
paquet de spin. 

On déduit des relations [2], [3] et [4] la valeur du 
coefficient d’amplication 9’ de la polarisation nucléaire 
dans le fluide 


== ] 
QE LS 


lo 


PO 1; 1 pe 


avec 


N° 11 
et 

TT CE) LA L Ÿ LA LA : 

1 sir fiv) G(v) dv Het 


Si l’on admet que la fonction de répartition G(w) a 
une forme de Lorentz et si l’on pose : 


GE 
le) LH En TAN — vo) 


on obtient l’expression suivante pour le coefficient 
d'amplification 6” dans le cas où y = vo + x 


ol Ta Vo ais(i + as}? Dis) (6) 
la VN ÿ 
avec 
1 3B [IT 
DIS 1e Fr (2) + as)? 
9 DE E1>\S B TES 
| PAT fl LE 7s)3/2. 
ee (1 Hbies (1 + a’s) 
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La formule (6), comme la formule (1), est valable 
dans le domaine ou s n’est pas trop élevé et où les 
termes d’ordre supérieur à s!/? dans ®’(s) ainsi que 
dans ®(s) sont négligeables. 

Dans le cas étudié ici, l’existence d’un temps 
d'accrochage vient encore limiter la rapidité de crois- 
sance de pb’ en fonction de l'intensité Z, du champ 
magnétique tournant. 


Lettre reçue le 12 août 1960. 
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LA RÉSONANCE | 
PARAMAGNÉTIQUE ÉLECTRONIQUE | 
D'UN MONOCRISTAL DE GLYCOCOLLE | 
IRRADIÉ PAR DES RAYONS GAMMA 


Par G. LomaGrio et J. UEBERSFELD, | 


Laboratoire de Spectroscopie Hertzienne, | 
Université de Besançon. 


Le monocristal étudié a été placé dans la cavité d’un 
spectromètre à résonance paramagnétique électro- 
nique de manière que l’axe de symétrie du monocristal 
(communément appelé l’axe b) coïncide avec l’axe de 


la cavité cylindrique et soit perpendiculaire au champ | 
magnétique constant /7,. | 


à 
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En collaboration avec E. Erb, l’un de nous a re- 


[marqué [1] que la forme du spectre d'absorption de ce 

J monocristal varie lorsqu'on le fait tourner autour de 
F 2e 

{son axe. D'autre part, les spectres accusent deux 


formes caractéristiques nettement différentes comme 


le montrent les photographies ci-dessous et il est diffi- 


Tete 


PIC: 


cile de penser qu’un seul radical hibre piégé dans le 
cristal soit responsable de ces deux types de spectre. 
Aussi avons-nous pensé [2], que la rupture de la liaison 
provoquée par l’irradiation donne naissance à deux 
radicaux libres, que nous avons tenté d'identifier. 


Essai d’identification des deux radicaux piégés dans 
le cristal. — Les différents spectres obtenus nous 
montrent l’existence de deux groupes de triplets laté- 
raux, symétriques par rapport à l’axe du spectre, qu’il 
est tentant a priori d'attribuer à un groupement NH ; 
mais nous pensons comme Ghosh et Whiffen [3], que 
le couplage entre le spin électronique et le noyau de 
l’atome d’azote et les trois protons est insuffisant pour 
rendre compte de l'importance de l’espace noté (bf) 
sur la figure 3. 

C’est pourquoi notre attention s’est portée sur cha- 
cune des mesures des intervalles (ab), (fg), (ag) à 
partir de tous les spectres obtenus. Cet examen nous a 
montré que dans tous les cas : 

nie le) — *: ou encore : 


(bi) =" 3 (ag) 5. 
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Autrement dit, nous pensons que les groupes des tri- 
plets latéraux (chaque triplet étant dû à l'interaction 
hyperfine avec l'atome d'azote) font partie d’un spectre 
de six raies équidistantes, que l’on peut attribuer à 


d 


IN, 


l'interaction d’un électron célibataire avec 5 protons 
également couplés. 

Nous avons donc supposé, qu’un des radicaux piégé 
dans le cristal donnait un tel spectre. (11 peut s’agir en 
effet d’un radical de la forme *NH, — CH, puisque la 
formule du glycocolle cristallin est 

NERO CENT 

la photographie figure 2, le laisse d’ailleurs pressentir. 
Nous avons reconstitué théoriquement le spectre de ce 
radical en donnant à chacune des raies la largeur 
moyenne repérée sur les photographies et les intensités 
convenables, qui pour notre hypothèse de 5 protons 
également couplés sont proportionnelles aux nombres 1, 
DD O 51e 


AH 


LHC 


Cette reconstitution a été réalisée pour chaque 
spectre obtenu, et la différence entre le spectre observé 
et le spectre reconstitué du radical NH, — CH, nous a 
fourni un ensemble de spectres, toujours en forme de 
triplet, évoluant rapidement en fonction de la rotation 
du eristal. La figure 4 illustre notre travail d'analyse, 


804 JOURNAL DE PHYSIQUE 


qui une fois terminé nous a permis de penser que les 
spectres observés sont formés de la superposition d’un 
spectre de six raies équidistantes pouvant être attribué 
au radical *NH, — CH, et d’un triplet attribuable à 
un radical CH ou CH, — CO;. 

A partir de cette analyse, nous avons construit les 
courbes représentant les variations des séparations AH 
et AH des raies de chacun des spectres en fonction 
de la rotation du cristal. Les variations sont sinu- 
soidales. 

Nous avons vérifié que l’on peut rendre compte de 
ces variations à partir de l’hamiltonien de spin de radi- 
caux libres piégés dans un cristal monoclinique. En 
partant de l’hamiltonien W (W = g8HS + F;), 
où g 5 HS est l’action directe du champ magnétique 
sur le spin et F, la fonction perturbatrice due à l’inter- 
action entre les spins électroniques et les noyaux, et 
en développant au premier ordre, nous avons obtenu 
l'expression suivante caractérisant les différentes fré- 
quences v des transitions entre les sous-niveaux 
Zeemann : 


NN = 50 FIRE fo) Z mi 


(m; représente le spin d’un proton et o l’angle dont 
tourne le cristal) ; la fonction f(o) s'écrit : 


f(®) — # cos 29 + B sin 20 + C}, 


Dans ces formules À désigne la constante de Planck, 
c la vitesse de la lumière, 8 le magnéton de Bohret g le 
facteur de Landé. 

La courbe représentant /(@) ayant même allure que 
les courbes représentant les variations des sépara- 
tions AZJ et AH” des raies de chacun des spectres, nous 
avons réalisé les deux identifications suivantes : 


AH = fl?) et AH = j{). 


Ce qui nous a permis de calculer pour chaque radical 
les constantes #, B, C. Ces constantes ne nous donnent 
qu’un ordre de grandeur des couplages hyperfins à 
cause de l'hypothèse très simplificatrice à la base même 
de ce travail, qui consiste à admettre que les protons 
sont également couplés dans chaque radicaliibre étudié. 

Nos résultats sont les suivants : pour le radical CH: 
pour lequel nous nous sommes donné $ — 0, ce qui 
revient à faire un déphasage arbitraire, nous avons 
trouvé : 


CO 1l0RICHAEE CID COS 
et pour le radical NH, — CH, : 
HP — A0 cmt, CN 10cm. 


Ces valeurs sont raisonnables, car toutes inférieures 
à l'interaction qui existe entre le spin de l’électron de 
Phydrogène et celui de son noyau. 


Lettre reçue le 22 août 1960. 
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AU SUJET DE L'HAMILTONIEN DE SPIN 
APPLICABLE A UN MONOCRISTAL 
APPARTENANT 
AU SYSTÈME ORTHORHOMBIQUE 


Par G. LoMAGLio, 


Laboratoire de Spectroscopie Hertzienne, 
Université de Besançon. 


Le but de cet exposé est de donner l’expression géné-| 
rale des niveaux d’énergie d’un centre paramagnétique, 
calculés au premier et au deuxième ordre. L'intérêt 
d’un tel calcul est de permettre l’interprétation des 
spectres de résonance paramagnétique des « radicaux 
libres » créés par irradiation dans des substances orga- 
niques à faible symétrie. Nous nous plaçons dans le 
cas où l'interaction spin-spin est anisotrope et où le 
noyau possède un moment quadrupolaire non négli- 
geable, mais nous négligeons l’action directe du champ 
sur le noyau. | 

Nous exprimons lPhamiltonien de spin en fonction 
des données suivantes : 

a) les trois angles d’Euler à, 0, ©, repérant la direc- 
tion du champ magnétique par rapport aux axes du 
champ cristallin (x, y, 2) ; 

b) les trois angles d’Euler uw, , v repérant l’axe de | 
quantification du spin effectif S'; 

c) les trois angles d’'Euler U, Ÿ, V repérant l’axe de | 
quantification des spins nucléaires. 

Pour un champ électrique cristallin à symétrie or- 
thorhombique, lhamiltonien de spin s’écrit [1] : 


Je — (gx HxSx + £y FR Sy + £z H,$;) 
Dis —5ss + 1) | + E(S2 — 655) 


PAS LEURS PISTES 


se P| PS II + | + P(TÈ— 1); (1) | 


où S et Z sont les opérateurs des spins électronique et 
nucléaire, H,, H,, H,les composantes du champ magné- 
tique dans le système de coordonnées x, V, Z; gx, £y, © 

les composantes du tenseur de susceptibilité, 6 le ma- 
gnéton de Bohr, D et Æ les constantes de structure 
fine, A,, B,, C, celles de structure hyperfine, P et P’ 

les constantes de couplage quadrupolaire. 

Le premier terme de J, est le terme principal, les 
termes suivants considérés comme des perturbations : 


1) 
D Sè—5 ss UN CRT (SANS) 
est la perturbation due au champ cristallin, les trois 
termes suivants constituent la perturbation due à 


l'interaction hyperfine, et les deux derniers sont relatifs 
à la perturbation due au moment quadrupolaire. 


1. Caleul des niveaux d’énergie au premier ordre, — 
a) Le terme principal de l’hamiltonien s’écrit plus sim- 
plement : 


Kp = 6 HogM (2) 


| 


| 


Il 


| 
| 


| 


g? — sin?0 sin°?3 gË + sin?20 cos?à g? + cos20 


ment d’axes : 


LS — 8 8H, M + M? 


2 
82) 
15 


ce qu’on écrit encore : 


g? = sin?0 gf + cos?0 gÿ 


2 : © 
gun — 82 et gi = sin? g£ + cos’? g?. 


b) Le calcul de la perturbation due au champ cris- 


| tallin nécessite la connaissance des éléments de matrice 
= MIS,M' >, < MIS! 
\ M' == 


ME < MISE 
M + Let l'application des formules 


M'> avec 
de change- 


C _… 4 7 
DE = y 9x + Bi Sy == Y19z 
a +7 <e LL 
Sy Lo 9x + Bay + Y2 923 


Sz = a3Sx + PsSy + YsSz 


Lafin d'exprimer la perturbation en fonction des quan- 
btités S’. 


La valeur moyenne du potentiel perturbateur dans 


| Pétat non perturbé nous donne le développement de la 
* perturbation au premier ordre : 


p re 
M2 F (3 cos? d — 1) ee cos 2u(1 — cos°? si] (4) 


pour le terme dépendant de MW?, le terme indépendant 
de M est laissé de côté, car 1l ne donne pas de contri- 
bution dans le calcul des fréquences des transitions. 


Notons encore que les angles d et u sont liés aux 
L angles 0 et Ô par les relations : 


: ‘à . 
SU — É sin 0 sin u — 2 sin D 
g 
5 5L 2 
(4 bis) 
L2Z L22 
cos L = 2 cos 0 cos u — ©? cos à. 
81 


L’hamiltonien de structure fine développé au pre- 
mier ordre s’écrit donc : 


(5) 
2 2 
gi cos? 0 
£ re | | 


c) Dans le champ de constant 5, chaque 
état de spin effectif marqué | M => consiste en 27 + 1 


de & 
— — COS 2u 
? 


D |3gjf cos? 0 1) 
2 g? 


} sous-états notés sous la forme |W, m >= au point de 


vue de la structure hyperfine (m# prenant les valeurs : 
PI — 1,7 —2,...,—1). 

_ Pour évaluer la perturbation due à l'interaction 
hyperfine, nous sommes amenés à évaluer les éléments 
de matrice : 


Tr ml$z LM, T2, 
CRT ms + I_|M +1,m—17, 


€ M, m|S + 1e [M +1,m+17, 


NS mis = 13 [M — 1,m+1>7, 
€ M,mlS= 1_|M—1,m—17, 


et à exprimer les produits : 
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A3 933 + Br Sx lx + CySy 14 en fonction des quan- 
tités connues S’/”. Ce travail effectué, nous trouvons 
que la perturbation de structure hyperfine développée 
au premier ordre est : 


Mm Ÿ A, cos Ÿ cos Y 
+ B, sin u sin U + Cy cos u cos U) sin 4 sin Ÿ L. (6) 
En faisant dans cette expression w DA: 
By; = Cy = B, À, — À, nous retrouvons le résultat 


classique relatif au champ électrique cristallin à symé- 
trie axiale : soit À Mm avecles formules de Bleaney [2]: 


K? 9? — A? gf, cos? 0 + B? gi sin? 0 
, Bg; sin 0 Ag, cos 0 
se ps DSL à up __ “18n COS Fe 
sin NS Te et cos PET RR (7) 
o o 


d) Le calcul de la perturbation due au moment qua- 
drupolaire est très simple, il se fait d’une manière 
analogue à celui de la perturbation due au champ 
cristallin. Il nécessite la connaissance des éléments de 
matrice : 


< mllzhn ee 
<< ml m 


Pour le développement au premier ordre de cette 
perturbation, le résultat obtenu est le suivant : 


ed + 


K°? g? 


< m|1 . im' > 


avec m'—=m +. 


2e 7 gi sin? 0 
F 2 2 
K° 8 


(8) 


et un terme indépendant de m», que nous laissons de 
côté. Remarquons immédiatement que ce terme est 
nul pour une transition de la forme 

(M, m—(M —1, m). 

Les paragraphes précédents c) et d) nous donnent, 
par les formules (6) et (8) compte tenu de (7), les 
termes à ajouter à l’hamiltonien de spin donné sous 
la forme (5) lorsqu'il y a interaction hyperfine ét un 
moment quadrupolaire, les perturbations n’étant déve- 
loppées qu’au premier ordre. 


2. Extension au deuxième ordre. — Nous avons 
poussé le développement de ces trois perturbations 
Jusqu'au second ordre, mais les calculs qui permettent 
d'aboutir aux résultats sont très longs et nous ne pou- 
vons les reproduire tous dans cette note. En particulier 
celui du développement au second ordre de l’inter- 
action hyperfine est une fonction compliquée de tous 
les paramètres u, Ü, v, V, 4, W', qui tient à elle seule 
plusieurs pages et qui sera publiée ailleurs. Quant aux 
développements des perturbations dues au champ cris- 
tallin et au moment quadrupolaire, ils sont plus 
simples et nous nous bornerons à les présenter par les 
formules suivantes (9) et (10) : 


1 
88H 


2M(2M? + 1[&) —2MS(S +1) 


gi 81 
| ILBL in? 6 cos? 0(D + E cos 2u)? 


6 
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TAUSITZ NA EEE 
“4 [£ =) E? sin? | 


(2 
o 


se MR t2MS(S + 1) —2M M 


sg PH, 


à cos ? 0 cf cos? 0) |” 
| [2(: ARR ) E cos (1 + Hoenl 
CE £ Ë 


: op COS 0\° 
L 4 (£ sin 2u 2e) | (9) 


0 
o 


Il 


4ami(l + 1) — 8m — m 
2KM(M — 1) 


AB gy g1 sin 0 cos JL 


De 57 9 I]\2 
(1 PECOS2U) ( Ke gt 
RP Sin nt gi sin? 0 } — m[21(1 +1) — 2m? — 1] 


K? 9? | SKM(M — 1) 


A? gjf cos? ‘)] È 


{ PB? gi sin? 0 
| KE 


| PPT P'o08 QU (1 = 


12 9 cos? 0 | 


+ 4P°?$in? 2U- (10) 


L'expression (9) est conforme à celle donnée par 


Weger et Low. L'expression (10) donnant l'influence 
du moment quadrupolaire n’a jamais été publiée à notre 
connaissance. Avant de terminer nous pouvons encore 
signaler que nous retrouvons tous les résultats établis 
par Bleanev à partir des nôtres, car pour nous placer 
dans le cas d’un champ cristallin à symétrie axiale, 
il suffit de faire dans nos formules : 4 — U — 0, 

NV 0 RE OT PES D ASS ANR = CB; 

Nous tenons à remercier ici M. Uebersfeld, Profes- 
seur à la Faculté des Sciences pour ses conseils et son 
aide. 

| Lettre reçue le 22 août 1960. 
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MESURES DE CORRÉLATION BÉTA-GAMMA 
POLARISATION CIRCULAIRE SUR l’SAu ET !2Sp, 


Par J. P. Deursex (*) et P. Lapwix (*) 


Centre de Physique Nucléaire, 
Université de Louvain, Belgique. 
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Les résultats préliminaires de ces mesures ont été 
communiqués au Colloque de Grenoble de la Société 
Française de Physique (29 février, 1er et 2 mars 


(*) Chercheur agréé de l’I. L. S. N., Belgique. 
(*#) Des calculs et des mesures d’analysibilité du polari- 
metre sont en cours et seront publiés ultérieurement. 
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1960) [4]. Nous avions signalé à cette occasion une In: 
fluence du champ magnétique de fuite du diffuseur sur] 
le circuit de coïncidence et nous avions évalué la cor: 
reetion à apporter aux résultats. 

Dans cette lettre, nous publions les résultats obtenus 
en éliminant cette erreur systématique. A défaut d’étas 
lonnage absolu (**) nous rapportons la polarisation 
trouvée dans 12?Sb à celle mesurée pour 1’Au ; nous| 
acceptons comme polarisation réelle la moyenne, pon:} 
dérée par leurs erreurs, des résultats rapportés park 


d’autres chercheurs. Les résultats sont repris dans lef 
tableau 1. 


TABLEAU 1 


EXPÉRIENCE A; 
F,. Bogam et A. H. WapPsrTRA [4] 0,52 + 0,09 
J. Berraier et al. [5] 0,34 + 0,05 
R, M. STEFFEN [6] 0,45 + 0,07 
Valeur moyenne pondérée 0,40 + 0,04 


Ce procédé revient à « étalonner » notre polarimètrel! 
par la mesure de la polarisation circulaire du rayon! 
gamma de 411 keV de 1S#Au. Pour tenir compte de la 
différence d'énergie entre ce rayon gamma et celui 
de 122Sb (564 keV), nous multiplions l’analysibilité par 
un facteur correctif de 1,13 obtenu à partir des figuresk 
de mérite dela référence [3]. Nous négligeons donc la 
variation avec l'énergie de la contribution de la diffus} 
sion multiple ainsi que celle de l’aimantation effectivef 
éventuellement différente due à une pénétration plus 
ou moins grande dans le diffuseur. D’après les résultatsk 
préliminaires des calculs d’analysibilité, l'influence def 
la variation de ces facteurs est d’un ordre de grandeur 
inférieur aux erreurs statistiques. 

Le dispositif expérimental est décrit ailleurs [1, 2, 3L{ 
Sa bonne marche a été contrôlée : 1) en suivant l4} 
variation, due à la décroissance de source, du nombreb 
des signaux et du rapport des coïncidences réelles à 
celui des fortuites ; 2) en calculant la dispersion des! 
résultats partiels ct des résultats journaliers ; cette dis- 
persion se situait entre 0,7 et 0,8 écart standard. | 

L'influence du champ magnétique de fuite du polari- 


‘ mètre a été déterminée en contrôlant le nombre des! 


signaux dans les canaux latéraux (voir réf. 1) et en! 


effectuant des mesures pour une valeur 0 — 960 (0 est 
l'angle moyen formé par les directions d'émission des 
rayons bêta et gamma détectés) ; l'effet moyen obtenu 
a été de (0,036 + 0,115) %,. | 

Dans les mesures de 1°$Au la source était une feuille! 
d’or métallique de 270 microgramme/em? ; dans celles! 
de **#Sb elle était une feuille d’antimoine métallique! 
de 4 mg/cem?. (L’épaisseur du support de source était| 
négligeable par rapport à celle de la source.) | 

L'influence d’une diffusion éventuelle des rayons} 
bêta dans les matériaux de Ia source a été contrôlée! 
de la façon suivante : nous avons couvert les deux côtés | 
d’une source de 1%#Au : 1) de deux feuilles d’or de 270 
microgrammes/cm? ; 2) d’une feuille d’antimoine de 4| 
mg/cm?. | 

Les résultats obtenus dans ces conditions sont, aux| 
erreurs statistiques près, identiques à ceux obtenus! 
avec les sources de 18Aun non recouvertes. (Les effets | 
rapportés à ceux obtenus par les sources non recou- 


De Ne 1 


De vertes sont : 1) 1,10 0,46 et 2) 0,89 + 0,18.) 
t Remarquons d’ailleurs qu'une dépolarisation éven- 
| tuelle devrait être encore plus petite pour les rayons 
| bêta du 2#%Sb d'énergie moyenne plus grande que ceux 
_ de Æ8Au employés aux mesures de contrôle, 
. La diffusion dans l’air a été trouvée négligeable : il 
n’en était pas de même de la rétro-diffusion sur le dia- 
phragme de plomb situé non loin derrière la source 
_ (voir par ex. figure 1 de la référence 3) : cette rétro- 
diffusion a été éliminée en couvrant le diaphragme d’un 
matériau léger. 
A, est liée à l’effet e mesuré expérimentalement par 


[ 


LETTRES À LA RÉDACTION 807 


la relation (pour les notations voir par ex. référence 1) 


€ (4 


A, = = — 
717 9fcos 6vle Dot —P)  Kvjel — P) 


Utilisant la valeur moyenne pondérée de À, indiquée 
au tableau {, et nos mesures de corrélation bêta- 
gamma polarisation circulaire, nous obtenons pour la 
constante d'étalonnage À — 2f la valeur moyenne de 
K = (3,98 + 0,44) %. 

Pour 12?2$Sh nous obtenons les valeurs consignées au 
tableau 2. 


TABLEAU 2 


e vlc 


En 0:021+ 0,12 


0,890 + 0,025 


Cette valeur de À, nous fournit dans l’approxima- 
tion & la valeur du rapport V/Y pour la transition 
. bêta principale de ce noyau (pour la notation employée 
voir par ex. référence [7]. Nous obtenons : 

VJY = + 0,25 Æ 0,05. 
Une seconde solution serait V/Y > + 12). 

Nous venons d’être informés par R. M. Steffen de 
mesures analogues entreprises à Purdue University. 
Nous le remercions de ses informations et des dis- 
cussions intéressantes que nous avons pu avoir avec 
lui. M. le Pr M. de Hemptinne a bien voulu nous 
® montrer l'intérêt qu'il portait à ce travail ; nous avons 
+ largement profité des conseils de M. le Pr P. C. Macq. 
Nous tenons à leur exprimer nos plus vifs remer- 
ciements. 

Lettre reçue le 21 septembre 1960. 
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NOUVEAUX TYPES D'’'INTERACTION 
EN THÉORIE NON LOCALE 
DE CHAMP MÉSIQUE 


Par E. ARNOUS, 
Institut Henri-Poincaré. 


Un certain nombre de faits expérimentaux suggèrent 
l'existence d’une « coupure » des interactions au voisi- 
nage de la masse A du nucléon : par exemple la diffé- 
rence de masse x+ — x° [1], la différence de masse 
p—n [2], la diffusion électron-nucléon 131, etc... Il 


0,960 + 0,005 


K (2$b) — 


Rey 1.13 9, (AU) 


A, 


4,50 + 0,50 


— 0,034 +. 0,034 


est possible d'introduire cette coupure dans une théorie 
non locale par l'intermédiaire d’un facteur de forme 
p(p, 49) CG — P —T et Lorentz — invariant. Il suffit 
de choisir pour p une fonction convenable de l’invariant 


Q(p, 4) = p° 9° — (p-4)”, 
par exemple 


e(p, 4) = 2 m°?[(X2m? + Q), 


où À est la coupure — M et m la masse de la parti- 
cule [4]. - 

Or, dès que l’on introduit un facteur de forme, le 
nombre d’invariants que l’on peut utiliser dans l’hamil- 
tonien d'interaction augmente. En effet, en dehors des 
termes usuels d'interaction, « coupés » au voisinage 
de M, tout autre interaction est à priort possible, si 
elle n’agit qu'aux énergies voisines de #7, où l’expé- 
rience suggère « qu'il ne passe quelque chose ». De 
telles interactions sont actuellement inobservables en 
électrodynamique, étant donnée la précision expéri- 
mentale (la coupure est trop haute). Mais, en théorie 
mésique, elles peuvent jouer un rôle important et sont 
susceptibles d'expliquer des phénomènes tels que les 
résonances secondaires dans la diffusion pion-nucléon. 

Les différents types possibles d'interaction sont les 
suivants : En électrodynamique 


Hit. > ] di p d# a 4% k % (p) Tulp. a) V4) Aulh) 


= 4(p, g) vu + ib(p, a) (p + gju + cs(p, g)(p + qhu v”, 
où a correspond au terme usuel (a = p). 
best un scalaire et c, un pseudo-scalaire et 


me ee 


$ met QNimi Es Q 


Le premier facteur «éteint » l’interaction aux éner- 
gies un peu supérieures à À, le second facteur, 
QIC2 m? + ©), aux énergies inférieures à À. 

En théorie mésique, 


TE fa p dé q dt k Ÿ (p) Tip, 4) (a) D (k), 


D, = ap, q) v° + bip, 4) (p + qju Y vu 
+ c;(p, 4) + ds(p, 9) (p + qlut vu 
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a (ou b) correspond au terme usuel avec facteur de 
forme : a (ou b) = p. 
Les autres facteurs de forme sont du type 


AM? QICE M? + QP ; 
ils s’annulent à la limite locale À — co. 
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MESURE CALORIMÉTRIQUE 
DE LA PÉRIODE DE 1%#Au 


Par J. ROBERT, 
Laboratoire Curie. 


Nous avons mesuré la période radioactive de 1%8Au 
avec notre installation de microcalorimétrie adiaba- 
tique à compensation automatique décrite précé- 
demment [1]. 

L'or 198 décroit essentiellement par émission 8 de 
0,960 MeV, suivie d’une émission y de 0,412 MeV. 
Deux autres transitions $ de faible pourcentage abou- 
tissent l’une à un état excité de 1,088 MeV, l’autre à 
l’état fondamental. 

Divers travaux ont été déjà publiés sur la période 
de 1%#Au ; les principales valeurs obtenues sont, par 
ordre chronologique, les suivantes : 


Not | 


Silver (1949) T =12,69 +004 7°" 

Sinclair et Holloway (1951) T = 92,73. + 0,0 j SM 

Locket et Thomas (1953) T — 2,697 + 0,003 j [4] 

Bell et Yaffe (1954) T = 2,699 + 0,003 j [5] 

Tobailem (1955) T — 2,686 + 0,005 j [6]. 

Johansson (1956) T = 2,697 + 0,005 j {A 
(1956) 


Sastre et Price T' = 2,694 + 0,006 jJ [8] 

Ces différents auteurs ont utilisé pour ces mesures 
des compteurs G. M. ou des chambres d’ionisation. 

Nous avons fait irradier une feuille d’or de 192 mg 
à la pile EL 3 et avons suivi sa décroissance pendant 
cinq jours dans un calorimètre en or de 3 mm d’épais- 
seur pesant 66 g, absorbant la totalité du rayon-. 
nement & et une part importante du rayonnement Y. 
L'intensité de la source, au début des mesures, était de 
l’ordre de 2 curies ; l'élévation de température horaire 
était de l’ordre de 1,9 degré. 

La source a été passée au photo-tube après une 
décroissance suffisante :aucune rayonnement étranger 
(2®%Au, etc.) n’y a été décelé de façon appréciable. 

L'analyse des points expérimentaux par la méthode 
des moindres carrés a donné, comme valeur de la : 
période : 

: T = 64,79 + 0,09 heures : 
SOI 

2,699 + 0,004 jours. 


Cette valeur est en accord avec les précédents résul 
tats bien qu’obtenue par une méthode toute différente- 
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